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In certain 1+1 dimensional field theoretic toy models, one can go all the way from 
microscopic quarks via the hadron spectrum to the properties of hot and dense 
baryonic matter in an essentially analytic way. This "miracle" is illustrated through 
case studies of two popular large TV models, the Gross-Neveu and the 't Hooft 
model — caricatures of the Nambu-Jona-Lasinio model and real QCD, respectively. 
The main emphasis will be on aspects related to spontaneous symmetry breaking 
(discrete or continuous chiral symmetry, translational invariance) and confinement. 
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1 Introduction 

With the shift in emphasis in field theory from perturbative to non-perturba- 
tive phenomena during the last decades, exactly soluble models have acquired 
a quite respectable status. The limited number of analytical tools available at 
strong coupling, so painfully felt in quantum chromodynamics (QCD), makes 
it tempting to look for other sources of inspiration. Except for the celebrated 
case of supersymmetric N ^ 2 Yang-Mills theoryd, "soluble model" in quantum 
field theory still means "1+1 dimensional model". This makes it clear at once 
that we should not expect such toy models to be in any sense realistic; they 
can teach us only little (if anything) of direct relevance to the real world. 
Nevertheless, we find it quite appropriate to include them in a book devoted to 
analytic approaches to QCD. One can perhaps picture the relationship between 
QCD and the soluble toy models as follows: If QCD is the sun, then the 1+1 
dimensional models are like planets which are gravitating around it, shining 
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in its light. Their distance to the sun may vary and each of them develops a 
life of its own, although powered by the^sun. Only then can one understand 
why a theory like the Schwinger modelH describing massless electrons which 
only move along a line, have a linear Coulomb potential and form non-interac- 
ting, massive bosons can accumulate close to 1000 citations in the hep-archiv£ 
up to now. This striking success is shared by the Gross-Neveu (GNl modelH 
(massless fermions with a point interaction) and the 't Hooft modelo (QCD2 
in the large N limit). Whether one likes these models or not — they have 
become an integral part of the intellectual endeavor to understand the strong 
interactions and their traces pervade the QCD literature. 

There is another aspect which we should like to mention. A soluble model 
in quantum mechanics is typically solved once and for all, then put into the 
textbooks where it ceases to play any scientific role. By contrast, a soluble field 
theory is a laboratory which can be used over and over again for new types of 
questions — this is the reason why it is so much richer. When browsing through 
the numerous citations of models like the GN model or the 't Hooft model, 
one finds the whole history of QCD mapped out - ranging from questions of 
asymptotic freedom, quark-hadron duality, spectroscopy, scattering, over issues 
of confinement, covariance, vacuum structure, chiral symmetry breaking, up 
to sum rules, supersymmetric extensions, "stringy" aspects, heavy quarks, and 
finally hot and dense baryonic matter, our main topic here. Also some technical 
developments like light-cone quantization would be unthinkable without 2d 
models as firm testing ground. Whether in the end these models will help to 
solve the real hard problems in QCD is still an open question; they certainly 
help to prepare us intellectually and widen our horizon about the fascinating 
world of non-perturbative physics and non-Abelian gauge theories. 

As will be seen later on, one of the main themes of this article is chiral 
symmetry and its spontaneous breakdown (SSB). We are aware that study- 
ing SSB in low dimensions is like dancing on thin ice — there is the constant 
danger of falling into the trap of some no-go theorem and making a fool of 
oneself. Without claiming mathematical rigor, we believe that the large N 
limit provides a reasonable protection against these perils. Thereby, it signifi- 
cantly enlarges the spectrum of questions which can be addressed within lower 
dimensional soluble models. 

In the spirit of what has been said so far, we have attempted to write this 
article in a rather self-contained and hopefully pedagogical way. The literature 
about two-dimensional field theory is not always easy to access for the non- 
expert, mainly because a number of different, highly specialized techniques are 
employed by different sub-communities (e.g. bosonization, light-cone quanti- 
zation, functional methods, conformal and topological field theory). We have 
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refrained from using genuine 2d-technology as much as possible and instead 
based our presentation upon a simple yet powerfuLtcchnique familiar to ev- 
ery physicist, the Hartree-Fock (HF) approximationErQ For the large N models 
which we will consider here, one can actually go quite far with a relativistic 
extension of this basic tool from many-body theory. This should also make the 
article more readable for those who want to learn field theory while working 
on QCD related problems (e.g. students or researchers with a nuclear physics 
background like ourselves). 

Finally we should apologize for the perhaps not always balanced choice of 
citations: It is strongly biased towards work which fits into our approach and 
may not do justice to other important activities in the field. Our article is 
not intended as a comprehensive review article of the usual type, but rather 
an attempt to convey some of the intellectual pleasure of working on QCD 
motivated, analytically soluble problems, using as example of topical interest 
hot and dense matter in two-dimensional large N models. This is presumably 
more adequate for a contribution to a Festschrift like the present one in honor 
of Boris loffe. 



2 Basics about 1+1 dimensions 

The transition from the real world to 1-1-1 dimensions has drastic consequences. 
These consequences are responsible for tremendous technical simplifications of 
quantum field theory, but also for many unrealistic features and even patholo- 
gies which have to be kept in mind. In the present section, we have collected 
the main kinematical and geometrical facts which are independent of the spe- 
cific dynamics, as well as some characteristic aspects of two-dimensional gauge 
theories. The selection of topics is guided by the applications which we have 
in mind in later sections. Therefore, we leave out a large number of impor- 
tant topics which are closely associated with low dimensions but less relevant 
for usl3i3 These include boson-fermion mappings, existence of anyons with ar- 
bitrary statistics, conformal invariance, sigma models, integrability questions 
as well as relations to matrix models, string theories, and gravity. We will 
also make no use of the technique of light-cone quantization which has become 
a powerful tool for dealing with two-dimensional model field theoriesp It is 
not clear to us whether it is of any help in the study of field theory at finite 
temperature and chemical potential. 
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2. 1 Kinematics 

We use the metric and antisymmetric tensors 

in 1+1 dimensions. Clearly, the Lorentz group S0(1,1) consists only of boosts 
— there are no rotations (and hence there is no spin). The inhomogeneous 
Lorentz group or Poincare group IS0(1,1) has three generators (boost K, mo- 
mentum P, Hamiltonian H) which satisfy the Poincare algebra 

[H,P]=Q, [K,P]=iH, [K,H]=iP. (2) 

States are characterized by the Casimir operator 

P^P" = , (3) 



the dispersion relation has the standard form E = ^/p^ + rn? . A left-over type 
of helicity can be defined via Lorentz-boosts 



\ _ ( cosh^ sinh^ \ f 
J \ sinh^ cosh^ y I a;^ 



(4) 



with rapidity 

^ = artanh v . (5) 
States are said to have helicity s if they transform according to 

\s)'=e^^\s). (6) 

Lorentz scalars, vectors and tensors can be constructed in the usual way. If 
y is a two- vector, its light-cone components transform irreducibly according 
to helicity ±1, 

The trace of a tensor Tj^ is Lorentz scalar, as is an antisymmetric tensor (e.g. 
the electric field strength). A symmetric, tracclcss tensor has helicity 2, 

(TOO ^ yoiy ^ e±2e(rOo ^ yoi) 

One peculiar feature of massless states in two dimensions is the following: Any 
number of massless non- interacting particles moving in the same direction with 
momenta pi behave kinematically like a single massless particle of momentum 
YliiPi- This is partly responsible for the proliferation of massless modes in 
certain models. 



2.2 Dirac fermions 

We start our survey of various aspects of two-dimensional field theory with 
Dirac fermions which will be the main actors in the following. Canonically, 
fermions are distinguished from bosons by equal-time anticommutation rela- 
tions, 

{Mx),^j{y)}=d,jd{x~y) . (9) 

Since there is no spin, a Dirac fermion has only two components (particle/an- 
tiparticle degrees of freedom). The free Dirac equation reads 

(i7^a^ - m)V' = (10) 

where the 7-matrices are 2x2 matrices satisfying 

{7^7'^} = 2g'^'' ■ (11) 
Unless stated otherwise, they will be chosen as 

The free Dirac Hamiltonian (in first quantization) is 

/i = 7^^ + 70771 (13) 

with 

7^=7V=( J _?) . (14) 

This representation is well suited for the chiral limit m — )■ 0. Free spinors 
(eigenvalues ±Ep with Ep = ^/rn^ + p^) are 



u{p) = 



cos6lp/2 \ ^ 1 / y/Ep+p 
sin6lp/2 J yJEp-p 



For later convenience, we have introduced the free Bogoliubov angle as 

6lp = arccot f — ) G [0,7r]. (16) 
Vto/ 

It appears if one diagonalizes h of Eq. ([l^) by a canonical transformation with 
the unitary operator 

(17) 



Up = exp |-y (aj^bp - h\ap)\^ 



for each mode where ap and bp are second quantized operators for right- and 
lefthanded fermions with momentum p. In the massless hmit, right- and left- 
handed fermions are projected out by 

1 ± 7^ 

Pr,l = . (18) 

Chiral transformations are defined as 

^ ^ e'"''V , (19) 

and represent a symmetry of the massless theory (note that 7^ can be defined 
in an odd number of space dimensions only). To see the meaning of right/left- 
handedness in 1+1 dimensions, we specialize the above spinors to to = 0, 
where 

9p = e{-p)7T (20) 

and hence 

-iP) = ( Z^'l ) (21) 



(0(p) denotes the Heaviside step function). Thus handedness in two dimen- 
sions is correlated with the direction of motion ("right- and left movers"). The 
Lorentz transformation is consistent with the above definition of kinematic 
helicity s — ±1/2 for right- and left-handed fermions, respectively, 

iIj{x) ^ e«T'/V(:r') , (22) 

with x' as given in Eq. (^) . Another important peculiarity of fermions in two 
dimensions concerns the vector and axial currents, 

j^^^Pjf^^P, j^ = ^7V^ (23) 

and their (partial) conservation laws, 

« = 0, d^,fX = ^^^i^^'i^ ■ (24) 

Due to the severe restrictions for 7-matrices in two dimensions, these currents 
have actually only two independent components, 

Jv — Ja ' Jv — Ja ' (25) 

or 

fA = ^""M^ ■ (26) 
This intimate relationship between vector and axial currents is the key for 
understanding the existence of massless mesons and baryons (see Sec. 3.3). 
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2.3 Scalar bosons 



Bosons are quantized canonically with equal-time commutation relations 

[T,[x)A{.y)]^-yKx-v) . (27) 

Whereas the free massive scalar field has no special features as compared to 
3+1 dimensions and satisfies the Klein-Gordon equation 

{d^d^ + m')(^^Q , (28) 

massless bosons are quite delicate. The simplest object in 3-1-1 dimensional 
quantum field theory, the massless scalar field, becomea_the most subtle one 
in 1+1 dimensions due to severe infrared divergenceslij As pointed out by 
Colemanli^ formally, the integral appearing in the free two-point function 

(0|<^(x)</.(0)|0) = j ^e"=-<5(fc2)e(fco) (29) 

is IR divergent, as can be seen by performing the fc°-integration. 

Since this is also the simplest exanmle of a Goldstone boson in 3+1 dimensions 
(broken symmetry: ^ + c)li3 one expects trouble with the Goldstone 
theorem in two dimensions. This is indeed what happens: There is a rigorous 
no-go theorem which forbids SSB of a continuous symmetry in two dimensions 
(Coleman theoremO). A similar phenomenon is known in statistical physics: 
there is no long range order in two-dimensional systems (at least with short 
range interactionaij). The strong infrared fluctuations destroy the order, for 
instance in a crystal, a magnet or a superfluid. There is a different type of 
long range order though which we discuss in the next subsection. 

2.4 Long range order 

The absence of long range order in 1+1 dimensional QFT or two-dimensional 
statistical physics could be a fatal blow to our investigation — what interesting 
phase structure can possibly be left? To understand the way out, let us first 
come back to the essence of two-dimensionality alluded to in the preceding sub- 
section: Long range order cannot be maintained over arbitrarily large distances 
due to strong fluctuations of the Goldstone mode; there is no SSB. A nice ex- 
ample for this effect are atoms in a plane which want to form a crystalEll The 
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amplitude of fluctuations of each atom around its equilibrium position grows 
logarithmically with L, if L is the extension of the crystal. The calculation 
involves a similar infrared divergence as in the 2-point function (p9|), 

(A = 1/a, inverse lattice spacing). In order to circumvent this no-go theorem 
one can go to the large N limit, invoking a diverging number of degrees of 
freedom at each point. This was demonstrated by WittenEj'Ej for the chiral 
GN model, a case very pertinent to the present study. After integrating out 
the fermions with the help of a complex scalar auxiliary field $ — $oe'^, one 
gets the following effective low-energy action for the Goldstone mode 6, 

Aff = §^id,ef . (32) 
Chiral symmetry breaking can be probed by the correlator 

{Mx)^Hy)) ^ {'^*i^My)) ■ (33) 

It approaches a constant I'I'oP for — y| oo if the symmetry is broken. 
Without SSB, one expects an exponential decay 

($*(a;)$(y)) -e-^^l"-«l . (34) 

Here, a different behavior is found, 

($*(x)$(zj)) - (e-'^^^^e'^^^)) ~ e-2-G(— y)/^ _ . (35) 

\x - yl^''" 

G is the free massless boson propagator, logarithmic in two dimensions. Expo- 
nentiation of this propagator yields a power law behavior for the propagator. 
Evidently, there is no SSB at any finite iV, in agreement with the no-go theo- 
rem, but there is a loophole at TV = 00 since 

($*(x)$(y)) ^ 1 + ^ In \x -y\ + 0{llN^) . (36) 

The correlator becomes more and more long-range until it cannot be distin- 
guished from a constant, in the limit jMj^i-oo. A similar phenomenon is known 
from two-dimensional xy spin-modelO'tS'Ej under the name of topological or- 
der. 
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2. 5 Gauge fields 

Gauge symmetry implies redundant variables; this is nowhere as clear as in 
a two dimensional world where, due to the absence of transverse directions, 
one might guess that all degrees of freedom of a gauge field are "pure gauge" . 
This is not quite true — there are left-over quantum mechanical degrees of 
freedom stemming from gauge invariant zero-mode fields. They can only be 
discussed reasonably well on a finite interval of length L. Evidently, pure gauge 
theories are a totally different story in 1+1 and 3+1 dimensions, as are the 
corresponding interactions between static charges. 

2.5.1 Abelian gauge fields 

For the U(l) gauge theory, Lagrangian and field strength tensor look as usual, 

C - -i^^.F'^" , F^, = df,A, - d,A^ . (37) 

In two dimensions however, an antisymmetric tensor has only one independent 
component, a Lorcntz scalar (the electric field); there is no magnetic field. 
Canonical quantization in the Weyl gauge (Ao = 0) reveals that the electric 
field E = OqA^ is conjugate to (minus) the vector potential A^, 

[E{x),A'{y)]^iS{x~y) . (38) 

The Gauss law in this gauge has to be implemented as a constraint on the 
physical states (this is possible since G{x) is conserved), 

G{x)\) ^diE{x)\) ^0 . (39) 

This leaves only the 0-mode of Ai, the spatially constant electric field, as 
physical variable. Switching to finite L, we have 



y / dxA^{x) , e = / dxE{ 
L Jo Jo 



x) (40) 



i/^-i-e^, [e,a]=i. (41) 



with 

H = 

2L 

Since the integer part of g being the electric charge of static sources, can 
also be gauged away by a periodic gauge transformation 

U = , (42) 

pure QED reduces to quantum mechanics of a particle on a circleB If we 
couple static charges to the gauge field, this changes Gauss's law into G = 
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dxE — ep. The "longitudinal" electric field (i.e., the a;-dependent part) is the 
Coulomb field of the static charge, which is linear in one space dimension. For 
comparison with Yang-Mills theory below, assume two charges on a circle 
of length L at points x, y; then the change in energy is 

AE = g^Kix - y) (43) 

with the periodic Coulomb potential 

= -y^ + \\x~y\-^{x~yf. (44) 

Pure QED is confining in 1+1 dimensions for rather trivial reasons. In the 
limit L ^ oo, the fact that the potential to order is linear also follows on 
purely dimensional grounds, since [g] = L^^. Finally, we note that one has to 
impose the residual Gauss law {Q\) — 0); only the charge sector survives on 
the circle. 

2.5.2 Non-Ahelian gauge fields 

Cnnsirjer [hp Lagrangian of Yang-Mills theory with a SU(N) gauge group on a 
circlel 



21,1 



^ = —.K^F"^'' ' (45) 



and 



F% = d,A: - d^Al - gr'^A^Al . (46) 



Whereas canonical quantization in the Weyl gauge is as straightforward as in 
the Abelian case, the resolution of Gauss's law becomes more intricate due to 
its non-linearity 

Cix) = {DEYix) ^ d^.E" - gT^'^A^^^E'' . (47) 



Generically, the covariant derivative has — 1 zero modes; the projections of 
E onto these are the physical (quantum mechanical) variables e^cJ They are 
non-hermitean due to the projection onto A^-dependent basis vectors. The cor- 
responding gauge invariant coordinates are the eigenphases of the path ordered 
integrals (untraced Polyakov loops) around the compact space direction, 

r = Pe^-fo '^^^^'^^^ = Ve'S'^^V'' , a = oFt^ . (48) 
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One finds 

[ef,a9]=Wp, (49) 

and the Hamiltonian 

This form makes explicit the fact that the variables are curvilinear coordi- 
nates — Ha is just the Laplacian on the SU(N) group manifold, the Jacobian 
J^{a) being the reduced Haar measure (Vandermonde determinant for unitary 
matrices) 

J {a) J^sin^ C^gL{au - %j) j . (51) 
i>j ^ ^ 

Convenient angular variables are 

Lpi ^ gLaa (52) 

with the constraint ipi — for SU(N). If one introduces the analogue of a 
"radial wavefunction" 

the Laplacian is reduced to cartesian form (except for a center of mass correc- 
tion in the SU(N) case, as opposed to U(N)), and a constant effective potential 
appears. 

The reduced wave function u(lp) now has to satisfy the boundary conditions 

u(^) — if ifi^ ipj . (55) 

The ground state is simply $0 = const., normalizable owing to the compactness 
of the variables, with Eq = 0. To get the excited states for any N, it is advan- 
tageous tOpj^interpret this problem in terms of free, non-relativistic fermions 
on a circlecB (the Pauli principle is enforced by Eq. (^sf)) with periodic (an- 
tiperiodic) boundary conditions for odd (even) N. Then the wavefunctions are 
Slater determinants of single particle wave functions 

^„(^) = -^e™^ (56) 
V zvr 
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with m integer for odd, half integer for even N. The energy is 

where the sums run over aU occupied orbits. By filhng the lowest N single 
particle states, one can easily ascertain that the ground state wavefunction is 
constant and that the ground state energy vanishes. The configuration space is 
determined by the fundamental domain, the smallest region bounded by zeros 
of the Jacobian. We can choose 



Lpi <IP2 < ■■■ < fN (58) 

and the center of mass constraint J^i ^1 = for SU{N) (this represents a more 
complete gauge fixing). Since the ordering is on a circle rather than a line, 
it still leaves a residual gauge freedom, namely cyclic permutations of color 
labels. This residual Zpj symmetry is expected on topological grounds due to 

Ui{SUiN)/ZN) = Zn- (59) 

This is the well-known center symmetry xif pure Yang Mills theory which plays 
an important role at finite temperaturecB By way of illustration, take SU(2): 
There is only one independent angular variable, ip — {ipi ~ ip2)/'2 € [0,7r]. 
Then, YM2 reduces to quantum mechanics of a particle in an infinite square 
well with 

and 

Unif) ^ \-smrup , J{if) = sin^ If , = ^-^{n^ - I) ■ (61) 

V TT 8 

The Z2 center symmetry is just parity. A detailed discussion of the SU(3) case 
with two independent variables can be found in Ref. 27. At very large N , 
since the ipi repel each other (cf. Eq. ^5|)) and live on a circle, they will tend 
to distribute themselves along the circle like a pearl necklace, with fluctua- 
tions which are 1/A^ suppressed in amplitude. This observation will become 
important later on for the 't Hooft model at finite temperature. 

The existence of zero mode gluons also complicates the interaction between 
static charges, at least for finite L. If one puts a static qq-pair on a circle, the 
Coulomb potential evidently involves (covariant derivatives) rather than 
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d ^ , so that the Coulomb interaction gets modified by the gauge 
The Coulomb potential can be transformed to coordinate spacer'^^^' 



remnants. 



for i ^ j, whereas it is identical to the Abelian one, Eq. (^), for i ^ j. Since 
the (pi are dynamical variables, one is not yet through but still has to solve 
a fairly complicated quantum mechanical (N — l)-body problem with a two- 
body potential defined by the Coulomb interaction. This reflects the fact that 
due to color spin dynamics, charges are never static in the non-Abelian case, 
even if they arfi_nailed down in space. A full analytic solution is available for 
the 811(2) casep3 where a dynamical quantum mechanical supersymmetry is at 
workc3 It explains why the potential is strictly linear (at all L, not just in the 
limit L — > oo), and why - for vanishing separation of the charges - all excited 
states are doubly degenerate whereas the ground state is non-degenerate and 
has strictly zero energy, see Fig. |l|. In the limit L ^ oo, all excited states move 
up to infinite energy, the gauge field dynamics get frozen, and the difference 
between QED2 and YM2 disappears. 

The discussion about the zero mode gauge fields may seem a bit pedantic: 
Why bother about finite L at all? The reason is of course that we are interested 
among other things in the finite temperature problem. The physics relevance 
of the Polyakov loop is beyond doubt there. 

Finally, the fact that the left-over gluon variables are eigenphases of path- 
ordered integrals has been used in the literature to emphasize the string picture: 
Closed strings that wrap around the compact direction in pure Yang-Mills 
theory, open strings with quarks at the ends if one includes matter fields in the 
fundamental representationtSEj'O We shall not exploit this picture here any 
further. 

2.6 Renormalization 

From the kinetic term in the Lagrangian in two dimensions, one infers that 
Dirac fermion, scalar boson and vector fields have dimensions 

[^,]-L-V2, = (63) 

In self-interacting fermion theories, a quartic self-interaction has dimensionless 
coupling constant and is therefore perturbatively renormalizable (GN model) 
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L 

Figure 1: Interaction energy of static quarks at distance d on a circle of length L. Dashed 
Une: QED2, cf. Eq. Q; solid lines: SU(2) YM2 (cf. Refs. 28, 30). 

unlike in four dimensions (NJL model). Self- interacting bosons can have inter- 
actions of arbitrary high order with couplings of positive mass dimension and 
are all super rcnormalizable (e.g. sine-Gordon model). The gauge coupling 
has dimension of mass, therefore gauge theories like the Schwinger model or 
QCD2 are also super-renormalizable. 

In the U(l) theory, an axial anomaly occurs, as is well known from the 
Schwinger model. In one-flavor SU(N) QCD2, there is no axial anomaly. 

Pure fermionic theories of GN type are asymptotically free; they share this 
property with QCD and are rather exceptional in this respect. The /3-function 
is negative (see Sec. ^.Ll] ). If the bare mass is put equal to zero (as in the 
original GN model), dimensional transmutation occurs like in pure Yang-Mills 
theory in four dimensions: A massless theory generates mass dynamically. All 
dimensionless ratios can be predicted without free parameter, only the physical 
fermion mass is needed as an inputu 

In the massive GN model, there are two distinct parameters of dimension 
mass (matching the two bare parameters Ng^ , rno): The physical fermion mass 
m and the ratio iriQ / {N g'^) . The first parameter is directly observable via the 
spectrum. The second parameter is also observable: In the chiral GN model, 
it determines the mass of the would-be Goldstone boson, cf. Sec. |3.3[ In the 
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non-chiral GN model, it can be measured via the qq scattering length which is 
proportional to Ng^ /mo. In the limit toq — > 0, the pion mass goes to zero and 
the scattering length in the non-chiral model diverges. In this way, the second 
parameter disappears and all observables can be expressed in terms of m only. 

QCD2 is trivially asymptotically free, being super-renormalizable. The 
perturbative (short distance) interaction is determined by dimensional consid- 
erations. 

2.7 Soluble models of interest for QCD 

The first generatkin of soluble models like the Thirring modelH or the massless 
Schwinger modefl from ±1960 were not yet rich enough for the questions we are 
interested in. Although they can be solved exactly and have certain interesting 
non-pcrturbative aspects, they essentially describe free theories with a trivial 
S'-matrix. If one tries to make them non-trivial by adding e.g. a mass term 
or flavor, one loses solvability. The Schwinger model has confinement, but 
provides no model for hadrons. 

From this point of view, 't Hooft's idea of the 1/N expansiona is very 
nice indeed. By invoking a large number of fermion species and using 1/iV as a 
small parameter, one finds highly non-trivial results already to leading order. 
Semiclassical methods become applicable and yield systematic, reliable results. 
All symmetries like Poincare invariance, chiral symmetry, gauge invariance are 
preserved. The phenommjiological arguments in favor of the large N limit for 
QCD (Zwcig rule ctcr^rr^rl) naturally play no role in two dimensions, where 
it is a purely theoretical tool. 

Models based on the large N expansion can be considered as second gener- 
ation field theory models, also historically. The two models we^re particularly 
interested in came up in 1974, the GNEI and the 't Hooft modelel The GN model 
generalizes the Thirring model to N flavors, while the 't Hooft model general- 
izes the Schwinger model to the non-Abelian SU(N) gauge group with a large 
number of colors. They acquire a rich physical content, as we hope to convince 
the reader in the remainder of this article. Here, we only mention the following 
basic facts: The GN models describe self-interacting massless fermions with a 
point-like four-fermion interaction. Depending on the particular variant, they 
can have either a discrete or a continuous chiral syjnmetry; in the latter case, 
the GN model may be thought of as the NJL modeS in two dimensions. Mass 
terms which violate chiral symmetry explictly can be included as well without 
loss of solvability (at large N). The 't Hooft model is QCD2 with SU(N) color 
in the limit N 00, with massless or massive quarks in the fundamental rep- 
resentation. As mentioned above, the model on the infinite line has no gluon 
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degrees of freedom, but is also a self-interacting fermion theory, here with a 
linear Coulomb potential characteristic for the two dimensional world (in an 
appropriate gauge). Unlike the GN models, this model shows confinement, 
which makes it particularly interesting from the point of view of QCD. The 
fact that also in real life (in four dimensions), the NJL model and QCD are used 
in parallel nowadays for certain questions of strong interaction physics makes 
this pair of models a good match and immediately suggests many questions 
to be asked. Throughout this work, we do not consider any embellishments 
or generalizations (e.g. flavor, more complicated interactions, different matter 
fields) of these models which have occasionally been invoked. The nice thing 
about the original models is just the contrast between the simplicity of the 
Lagrangian and the complexity of the phenomena it produces, and we do not 
want to spoil this property. 

Finally, we should point out that the main goal of the present study are 
questions of hot and dense matter in these toy models. For this purpose, the 
large A'' limit is not only a technical trick leading to solvability of the models, 
but is in fact instrumental for enabling chiral symmetry breakdown in two 
dimensions. This is the reascp^sihy, another aspect of the GN model will play 
no role here, its integrabilityc3c3'El A lot of progress has been made towards 
solving the finite N model exactly. Here, we shall restrict ourselves to the 
large N limit so that the corresponding difference between the non-integrable 
't Hooft model and the integrable GN models will not play any role, but a 
common approach for both models will be used throughout. 

3 Two dimensional models at zero temperature and chemical po- 
tential 

3.1 Gross-Neveu models 

We cannot discuss finite temperature or density without first reviewing the 
particle content of the theory. Originally, theXJN model has been solved with 
the help of semi-classical functional methods^ supplemented by inverse scat- 
tering techniques from soliton theory for the baryonsEJ This approach can be 
rigorously justified in the large N limit. Typically, one introduces auxiliary 
bosonic fields and integrates out the fermions (Gaussian integral over Grass- 
mann variables), then applies a saddle point approximation to the remaining 
bosonic functional integral. More recently, an alternative to the invers e sc atter- 
ing method has been proposed, based on the Gel'fand-Dikii equationCJ Here, 
in order to have a coherent scheme for all questions that interest us as well 
as for pedagogical reasons, we shall make use of time honored concepts well 
known from atomic and nuclear physics, such as the HF and random phase 
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approximation (RPA), although in a fully relativistic setting0Ei0 In contrast 
to the functional integral methods, we thereby handle fermions more directly 
and, since we work canonically, can also more readily address questions of the 
vacuum structure (Dirac sea). Needless to say, the equations one eventually 
solves are always the same (HF or Schwinger-Dyson, RPA or Bethe-Salpeter 
equations). More than anything else the choice is a question of language. 

3.1.1 Vacuum, physical J ermion 

The Lagrangian of the GN model family, including a bare fermion mass term 
and suppressing color labels, readsEl 

C = qird^^q + [{qqf - A {qi^'qf] - m^qq . (64) 

The original model with discrete chiral symmetry q — > 7^9 is recovered for 
A = 0,TOo — 0; the choice A = l,mo = corresponds to the NJL-type model 
with continuous chiral symmetry; the toq term breaks both chiral symmetries 
explicitly. To leading order in 1/iV, the Hartree approximation can be used, 
replacing qq in the Euler-Lagrange equation 

17^9^9 + 5' {m - K(ll''Q)f) q-moq = (65) 

by its vacuum expectation value (the Fock term is 1 /N suppressed) . The non- 
zero value of the order parameter for chiral symmetry signals SSB. We first 
assume (qj^q) — 0, 

{i-f^df,-mo + Ng^Ps}q{x)=0 . (66) 
A physical fermion mass appears, 

m = mo~ Ng'^Ps , ps = ^(qq) , (67) 

as dictated by covariance. The essence of the Hartree approximation is to 
impose a self-consistency condition, here on the scalar density ps of the vacuum. 
Using the free fermion spinors of Eq. (|l5|) above we get 

(m - mo) , dfc 
Ng'^ 27r 

f+^/^ dk m m , A2 

= -/ — c.-— log— . 68 

J -A/2 27r + m2 27r 
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As usual, this is also a variational solution obtained by minimizing the vacuum 
energy in the space of Slater determinants. The Hartree- or gap equation, 



m = mo + ni—^ — In — ^ , (69) 

contains much information. For mp = (discrete or continuous model), it has 
two solutions, either m = or 

The broken symmetry solution is always lower in energy density by 

(--vac (--vac ~ ■ \''-) 

The gap equation for the broken phase does not determine m, as is clear from 
the fact that the theory contains no dimensionful parameter. If m is adjusted 
(to the "observed" fermion mass, whatever that means in two dimensions), 
then the gap equation is a renormalization condition which tells us how Ng^ 
(the bare coupling constant) depends on the cutoff A. Correspondingly, the 
physical mass has the dependence 

m = Aexp|-^} (72) 

familiar from similar dependences in real QCD and in line with the renormal- 
ization group. All physical quantities with dimension of mass have to depend 
in the same way on A. The essential singularity in 5 at 17 = underlines 
the power of the expansion as compared to ordinary perturbation the- 
ory in g. Dimensionless quantities can be predicted without any parameter, 
at mo = 0. Apart from this phenomenon of dimensional transmutationj— the 
other remarkable feature which we find here is asymptotic freedom (AF)eJ'Ea 
One can determine the /3-function of the GN model in an elementary way as 
follows: The physical fermion mass m should be cut-off independent, therefore 

AJL™(A,»(A))^(A^ + A^i|)„.0. (73) 



Set 



= m (74) 
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and find, by inserting m from Eq. (f3 



/3(5) = -^<0 (75) 

or AF (to tliis order, the result is independent of the renormalizalion scheme). 
More elaborate computations of the higher loop beta-functionsoo do not play 
a role in the large N limit. We shall return to the issue of AF in connection 
with qq scattering, where a running coupling constant will be defined. 

For mo = and the case of continuous chiral symmetry, the assumption 
(97^9) = is too restrictive; one gets a continuum of vacua related by U(l) 
transformations and labeled by a chiral angle (p, 

{m) = -^cosLp, (917^9) = ^ sin . (76) 
9 9 

For non- vanishing bare mass mg, the gap equation ( |69| ) signals the appearance 
of a second physical parameter -ff^, in addition to the overall scale set by m. 
Its relevance for scattering and bound state observables will become clear later 
on. 

3.1.2 Baryons 

It would be premature to conclude from the preceding section that the GN 
model at, iV — > cx) reduces to a free, massive fermion theory. As observed by 
Wittenp^l baryons have to be considered at large N and can be treated in HF 
approximation, like the vacuum. They correspond to a different kind of HF 
solution which breaks translational invariance. Their mass diverges like N, the 
baryon-baryon interaction also scales like N, and these effects have to be taken 
into account in leading order. Baryons weije-.found originally with functional 
techniques and inverse scattering methodsua Since we now know the scalar 
potential, we can use it to verify Witten's picture_pf baryons as relativistic 
HF solutions — this is what we propose to do herecj We consider the discrete 
chiral model first, since the continuous model has in addition totally different. 



light baryons, see Sec. ^ below. 

If we restrict ourselves to baryons made up of n < quarks, it is sufficient 
to take into account one positive energy "valence" level filled with n fermions, 
together with the completely filled Dirac sea. The Hartree equation without 
invoking translational invariance becomes 

{ird^.+Ng^p,{x)}q{x)^0 (77) 

where the scalar density now refers to the baryon state \B) , 

1 ^ 

1=1 
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We introduce the single particle energies and eigenfunctions of the Dirac Hamil- 
tonian corresponding to Eq. (|77|), 

7' Y ^ - I'Ng^p, {x) I Vi^) (x) = ±ea Vi^) {x) . (79) 

Then Ps{x) is given by 

Ps{x) = +Y.'^':^\x)i^i-\x) (80) 

Q 

where V'o^'' is the valence state. Dashen et azEl have found two different types 
of solutions for the scalar density, the "kink" and the "double kink" (or rather 
kink-antikink) . Let us start with the less exotic double kink, characterized by 
the following scalar potential: 

S{x) = -Ng^ps{x) = m{l + y (tanh^_ - tanh^+)} (81) 

with the definitions 

77- 77 

y^sinf?, (82) 

C± = ymx ± — artanhy . (83) 

S{x) has the form of a potential well and approaches for x ±00 the asymp- 
totic value m, i.e., the physical fcrmion mass. Varying y or equivalently the 
degree of occupation n/N, S{x) changes from a 1/ cosh^ (yma;) shape at small 
y to a Woods-Saxon like shape at j/ ~ 1, 

S{x) - (l + e(l"l-^)/^)"' , (84) 

with skin thickness c and radius R given by 

1 1 cos 6* , , , , 

2m 2to ^ 2 ' ' V ; 

respectively. We proceed to verify the self-consistency of this scalar potential. 
After a straightforward solution of the Dirac equation (^^ with potential (^l|), 
one finds the following results: The normalized wavefunctions of the discrete 
(valence) states are given by (in a basis where 7^^ = — cri,7^ = ias) 

4^\x) = ^ ( ) (86) 
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Figure 2: Single particle spectrum belonging to the baryon in the GN model. Continuum 
states fill the outer region \E/m\ > 1; dashed curves inside the gap are valence levels; y 
measures the occupation number of the positive energy valence state, cf. Eq. (63). 



with eigenvalues 



The continuum states are 



(±) 



= ±m\/l — 2/^ 



(87) 



4^' 



{x) = 



1 



V2E{k){ik + ym) 



(ifc — TO)(ifc 
±E{k)lik- 



ymtanh^-) 
ym tanh^-|_) 



with 



E{k) = y/k^ + m2 . (89) 
T}ipy are reflectionless, a characteristic feature of time independent solutions0 



U531 rpj^g gpectrum of the Dirac Hamiltonian therefore consists of the usual 
positive and negative energy continua starting from ±m and a pair of discrete 



states inside the mass gap at energies given by Eq. (j^), see Fig. g. These 
latter states are of course a prerequisite for obtaining bound baryons. The 
discrete negative energy state has to be filled completely {N fermions), its 
positive energy partner with n < N fermions (remember that, unlike in QCD, 
there is no restriction to color sing lets). Reinserting Eqs. (||,|8|) into (|l|), one 
can then verify the self-consistency of the solution. j-This requires the use of 
the vacuum gap equation (B9) to eliminate infinitiesua 
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Finally, one has to compute the mass of the baryon. This is somewhat 
delicate because of the vacuum subtraction. The single particle energies of two 
systems must be subtracted: One of the systems has only continuum states, 
whereas the other system has one extra pair of discrete states. The easiest way 
of doing this is to enclose the system in a finite box, thereby discretizing all 
states, and to use the well-known relation between phase shifts and density of 
states. In this way, one recovers the result of Dashen et al£3 for the baryon 
mass, 

2Nym / sme\ 

Mb - =nmi — \ . (90) 

Let us briefly compare this exact result to the non-relativistic approach 
(in the present model, the number of valence quarks governs the degree of 
relativity). In the non-relativistic Hartree approximation, we start from the 
Schrodinger equation for a particle of mass m, interacting with the other par- 
ticles through a (5 -function potentiaL. The Hartree equation then reduces to 
the non-linear Schrodinger equationpil 

^ ng^\cj,o\A M^) = EoM^) , (91) 



2rn dx'^ 

with the following bound state solution. 



1 

K^-nc/m, Eo^-—. (93) 
2 2m 



The baryon mass turns out to have the value 



Mb =nmll-l f ^ ) I . (94) 



6 V 2 
If we choose 

9^-^, (95) 

both Mb and the scalar density in the full calculation and in the non-relativistic 
limit agree up to 0{y^). The reason why the GN model corresponds to a 
specific choice of the coupling constant is that it involves only one parame- 
ter, whereas the non-relativistic model has two independent parameters, the 
fermion mass and the coupling constant. 
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Figure 3: Scalar potential S (solid), valence fermion density p'j (dashed) and valence energy 
Eg (dash-dotted lines) for GN model baryon and three different occupation fractions, n/N = 
0.5 (left), 0.9 (middle), and 0.99999 (right). In units of m. 



Now turn to the opposite, extreme relativistic limit. As approaches 7r/2, 
the effects of the Dirac sea become overwhelming. Eventually, the baryon goes 
over into a well separated kink anti-kink pair. Although the scalar density 
becomes flat, the fermion density for the discrete state is concentrated near 
the surface, since 

In the limit y — > 1, the two kinks are completely decoupled, see Fig. |[ 

If one looks at a single kink now, one finds another type of baryon of the 
GN model, referred to as Callan-Coleman-Gross-Zee kink in the literatureE3 
It has the scalar potential 

S'(a;) = m tanh (a;m) . (97) 

This solution has many unusual features: S{x) interpolates between the two 
different vacua related by the discrete 7^-transformation. Unlike the double 
kink baryon which is only stabilized by the binding of fermions, the single kink 
is topologically stable. The fact that a kink type potential can bind particles is 
a well-known, purely relativistic effect which has no analogue in non-relativistic 
quantum mechanics. The discrete state corresponding to Eq. (pq) becomes 



M^)^ ) , (98) 



with eigenvalue Eq = 0, while the continuum wavcfunctions replacing Eq. 
are 

^ ' ~ V2E{k) [ ±E{k) J ■ ^^^^ 
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The valence state has evidently vanishing scalar density so that there is no 
feedback at all to the Hartree potential. This is consistent with the fact that 
the mass of the kink baryon is independent of the number n of valence fermions. 
The self-consistency for the scalar potential can again be verified using the 
vacuum Hartree condition (|69|). As expected intuitively, the mass is 1/2 of the 
double kink mass in the limit y — > 1 

^Bkink - — , (100) 
TT 

in agreement with Dashen et a/Jll 
3.1.3 Mesons 

In this section, we shall rederive the known fermion-antifermion bound state 
(meson) of the GN model, using the relativistic RPA. This meson has zero 
binding energy in the large N limit, i.e., its mass is just twice the physical 
fermion massH We discuss only the discrete chiral symmetry case here, with- 
out bare mass; the NJL type model with its massless Goldstone mode will be 



approachOLj and solve them analytically for mesons with arbitrary momen- 
tum, thereby demonstrating explicitly the covariance of the spectrumo 

In RPA, the meson is described in terms of particle-hole excitations on 
top of the HF vacuum. We therefore start from the equation of motion of the 
following operator bilinear in the fermion fields, 

Qal3{x,y) = JJ'^ql|3{y)q^a{x) . (101) 
i 

Using the equation of motion for q, Eq. (|65|), and neglecting 1/N suppressed 
terms, we find 

i^g(x,2/) = -ii^^Q{x,y)^^ + ^^^Q{x,y)^ (102) 

-Ng^ {tr (7°Q(x, x)) 7°Q(x, y) - Q{x, y^hv {j°Q{y, y)) } 

Furthermore, we expand the operator Q{x,y) around its vacuum expectation 
value 

Q(x, y) = (0|Q(x, y)\0) + -^Q{x, y) + ... (103) 

V N 



treated iiJ-Sec. 3.3. We derive the RPA equations using an equations of motion 



where the c-number part corresponds to the density matrix of the vacuum, 

(0|Q(a;,y)|0) =p(a;-y) . (104) 
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Inserting (103) into (102) and equating terms with the same power in l/yiV, 
we get to zeroth order the following alternative formulation of the Hartree 
equation, 



i^p(-) 



4^[7^p(a;)] - 7VgV(7V(0)) b",p{x)] = (105) 



d ~ 



and to order l/vN the linearized equation of motion for the fluctuation Q, 

(106) 

-Ng^ {tr (7^(0)) j"Q{x, y) - Q{x, y)7°tr (7"p(0)) 

+ tr (7°Q(x, x)) 70p(x - y) - p{x - y)7°tr {j'^Qiy, y)) } . 

Switching to momentum variables, 

p(fc) = J dxe-'''''p{x) 
Q{k',k) = f dx [ dye-"''''+"'yQ{x,y) , (107) 



(108) 



the Hartree equation simplifies to 

[hik),p{k)] = [kr' +m-f°,p{k)] = 



where h{k) denotes the single particle Dirac Hamiltonian and where we have 
again introduced the physical fermion mass, Eq. (|67|). The Hartree solution of 
Sec. p.l.l| corresponds to 



p{k) = v{k)vHk) = - - 



2 2E{k) ' 2E{k) 



(109) 



and obviously satisfies Eq. (108). In order to convert Eq. (106) into the canoni- 
cal form of the RPA, we first have to project the operator Q onto self-consistent 
Hartree spinors. In the present case this simply amounts to a change of ba- 
sis from massless to massive free spinors. We define the "particle-hole" and 
"hole-particle" components of the operator Q as 



]i2(k',k) 
hi{k',k) 



u\k')Q{k' ,k)v{k) 
v\k')Q{k' ,k)u{k) 



(110) 



26 



Using the commutation relations, the corresponding "particle-particle" and 
"hole-hole" operators Qu and Q22 can be shown to be of higher order in 1/VlV 
and can therefore be neglected. The resulting system of coupled equations is 
sandwiched between the vacuum |0) and a one meson state {n,P\, where P 
denotes the total momentum and n the type of meson. We define the transition 
matrix-elements 

(n,P|Q2i(fc',A:)|0) = 2Tr5{P - k + k')Xn{P, k) , 
(n,P\Qi2ik',k)\0) ^ 2nS{P-k + k')Yn{P,k) . (Ill) 

Then, using the Heisenberg equation of motion in order to replace time deriva- 
tives by energies, we obtain 

£niP)XniP,k) = E{k~ P,k)Xn{P,k) (112) 

/Ay 
-^{v{k')u{k' -P)Y^{P,k') 

+uik')v{k' ~P)X„{P,k')} 
£n{P)YniP,k) = -E{k-P,k)Yn{P,k) (113) 

/dk' 
{v{k')u{k' - P)Yn{P, k') 
27r 



+u{k')v{k' - P)Xr,{P,k')} 



with 



E{k',k) ^ E{k') + E{k) . (114) 
These are the final RPA equations. Taking into account the symmetry relation 

u{k)v{k') = -v{k)u{k') (115) 

which holds for the spinors ( p^ ) , they can be cast into the familiar matrix form 
of the RPA, 

A B \ f X \ ^ f X 



-B -AJ^Y )-'{y) ■ (116) 

It is illustrated graphically in Fig. IJ. Characteristic of the RPA (as opposed 
to the Tamm-Dancoff approximation (TDA)) are the backward going bubbles 
(corresponding to y-components) which account for ground state correlations. 
Since the integral kernels are separable, one_can solve these equations alge- 
braically. After a few straightforward stepaij one arrives at the eigenvalue 
equation 

_Ng'^ /■ dfc / 1 1 \Am^ + P^~E^{k~P,k) 



2n\E{k-P) E{k)J £1{P) - E^{k- P,k) 
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Figure 4: Diagrammatic illustration of HF (upper part) and RPA (lower part) for the GN 
model. Dashed lines: bare (massless) quark propagator, solid lines: dressed (massive) quark 
propagator. 



Clearly, the choice 



(118) 



converts this into the vacuum Hartree equation. This shows that there is a 
marginally bound meson with mass 2m, and that the spectrum is covarianL 
It is instructive to compute also the qq scattering matrix in this approachlEj. 
The scattering matrix can easily be found, since we are dealing with a separable 
potential. With the Mandelstam variable s = 2P£, find 



V2 



-is) 



V^-V log 



i-yi^ 

l + VT^ 



(119) 



in the scattering regime with rj = Am? j s < 1. (r is normalized such that r = 
—Ng^m^/y/S in tree approximation). r(s) has the expected pole at s = 4rn^ 
corresponding to the threshold bound state. By going to large s, we can check 
asymptotic freedom in a very physical manner. 



t{s) 



1 



^ log(m2/s) 



(120) 



The interaction shows the logarithmic decrease in the UV characteristic for 
asymptotically free theories. A running coupling constant can be defined by 



28 



C = qilpq - -trF^.Ff^" , ^{9^ + igA^) . (122) 



comparison with the tree approximation (s = ^^), 

iV5cfr(M) = 1 2 ■ (121) 

This iUustrates once again the usefulness of the GN model as toy model for 
QCD. 

3.2 't Hooft model 

The 't Hooft modell is defined as the large N limit of 1+1 dimensional SU(A'^) 
gauge theory with quarks in the fundamental representation, 

1 

-1 

2 

In the axial gauge, the gluons are gauged away, leaving behind a theory of 
fermions interacting via a Coulomb potential which is linearly rising in one 
space dimension, 't Hooft originally solved this model for the meson spec- 
trum, using light-cone coordinates and the light-cone gauge. This was the firjat 
practical demonstration of the power of the 1 /N expansion for gauge theoriesL3 
Working diagrammatically, he identified planar diagrams without quark loops 
as being of leading order in 1/A^ (nested rainbow graphs and ladders) and 
summed them up by solving (light-cone) Schwinger-Dyson and Bethe-Salpeter 
equations. He thus found an infinite tower of mesons lying on a kind of Regge 
trajectory. Early applications of this model dealt. wMh 'if^tlPf of confinement, 
asymptotic freedom, form factors, scattering etcElOSEElO For some time, 
the precise character of the 't Hooft limit was somewhat disputed (strong vs. 
weak coupling, order of limits, 't Hooft boundary conditiao. IR regularization 
of the quark self-energy, Wu's alternative meson equatiorO). Studies in ordi- 
nary coordinates fnllnwpH wViprp thp vacuum and baryons can be dealt with 
in a more direct wayj£300Eje they fully confirmed the original light-cone 
results for the meson spectrum. During the last few years, various extensions 
of the 't Hooft model were fffupd useful as laboratory for QCD related ques- 
tions, such as adjoint quarksEEfEZl (a toy model fpr^uDns and strings), heavy 
quarks, operator product expansion, duality etcEjOu In line with our main 



subject, we propose to go through the same three topics as in Sec. 3.1 for the 



GN model — the vacuum, baryons and mesons. The light sector (mesons and 



baryons which become massless in the chiral limit) will be deferred to Sec. 3.3 
where we will cover this aspect in a unified way for both models. 
3.2.1 Vacuum, confinement 

We start with a brief reminder of the vacuum, cjhiral symmetry breaking and 
the role of confinement in the 't Hooft modeluflca Since this latter aspect is 
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intimately related to infrared singularities (absent in the GN model), we work 
in a finite interval of length L, taking the limit L — > cxd at the end whenever 
possible. Fermions are assumed to obey antiperiodic boundary conditions. We 
introduce a bilinear operator which describes color singlet dynamics to leading 
order in large N (cf. GN model in Sec. 3.1.3| ), 

The ai^hi are annihilation operators for right- and left-handed quarks. The 
Hamiltonian for the 't Hooft model in the axial gauge has the form 

^ = Y.T^P+^l'^){''\^P>^^P)-^l^P)^^^P)) (124) 



167r2 

The currents jij{n) can be taken in the U(A^) form at large N , 

j^J{n)=Y,(a]{p)a,{p + n) + h]{p)h{p + n)) . (125) 
p 

H is expressed in terms of Q as 

^ = ^^E (P+^)*^^'Q(^''P)+™^E*^^°Q(^''P) (126) 

N'^g'^L 1 
+ ^Q^2 E ■^ir{6rsQ{r,r)-Q{s,r)Q{r-n,s-n)} , 

whereas the equation of motion for Q becomes 

idtQ{p',p) = ^{{p' + l/2)Y'Q{p',p)^{p+l/2)Q{p',ph') (127) 

+m[7",Q(p',p)]-|-^ E ^{Q{p',P")Q{p"-n,p-n) 

-Qip'^,2,p"-n)Qip'\p)} . 
Expansion of Q around its c-number part (vacuum expectation value) 

Qip',p) - {0\Q{p',p)\0) + ^Q{p',p) + ... 

V A* 



Sp,ppip) + ^Qip',p) + ... (128) 
V J\ 
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Figure 5; Same as Fig. ^ but for 't Hooft model in axial gauge. Curly lines: Static Coulomb 
potential. The extra diagram in 't Hooft's double line representation (upper right) explains 
the dominance of the Fock term at large A'^. 



yields, to leading order, the HF equation 

[hHF,p{p)] = (129) 
with the single particle HF Hamiltonian 

hHF - ^(P + 1/2)7' + mj" ^ - ") - 1/2) . (130) 

n=iO 

In contrast to the GN model where the Hartree term dominates, here only the 
Fock term survives at large N, see Fig. || (the Coulomb interaction acts like a 
color exchange force). Solving the HF equation self-consistently is equijszalent 
to summing up nested rainbow graphs in a Schwinger-Dyson equationaH On 
general grounds, the vacuum density matrix p can be parametrized as 

P{P) = ^ +7Vo(p) -i7Vi(p) +7V5(p) (131) 

where pi is only allowed if m = 0. Charge conjugation requires po,i to be 
even, p^ to be odd functions of momenta. The Slater determinant condition 
p^{p) — p{p) characteristic for HF then translates into 

pI{p)+pI{p)+pI{p)^\. (132) 
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Consider first the case of massive quarks. For m 7^ 0, pi vanishes and we can 
set ^ ^ 

PQ = --sme{p), p5 = --cose{p) . (133) 

Comparing 

p{p) ^ v{p)vHp) (134) 

with the HF spinors of Sec. ^.2| (Eq. |l5|), 9{p) is recognized as the Bogohubov 
angle. The HF equation now becomes 

'^{p + 1/2) sme{p) ~ mcoseip) + ^^^^y^ ^ smie{p) - e{p - n)) = , 

n=iQ 

(135) 

i.e., a nonlinear integral equation in the limit L — s- 00. For m — 0, pi can 
be non-zero. As in the chiral GN model, there is now a continuum of vacua 
parametrized by a chiral angle if, 




^ I s\n9{p) coBip 
'n9{p) sin 
cos6'(p) 



- \ sin 6l(p) sin ^ | . (136) 



The gap equation (135) remains valid provided we set m — 0. The chiral 
condensates are 



Although thij^afl equation (|135| ) for the 't Hooft model has to be solved 
numericallyE^fEjO the value of the quark condensate is known analytically, 
owing to-an indirect determination via sum rules and the 't Hooft equation for 
mesons;0 it is (for — 0) 

(This formula has been generalized to the case of finite bare quark mass E) 
In Fig. ^, we show a few typical results for the Bogoliubov angle for several 
values of the quark mass and compare them to the free theory. For larger quark 
masses, the interaction effects are very small and can be treated perturbatively. 

So far, there seems to be little difference between the GN and 't Hooft 
models — superficially, the vacuum looks in both cases as if the fermions had 
acquired a dynamical mass, although momentum dependent in the 't Hooft 
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case. How does confinement of quarks manifest itself? In QCD2, there are 
no massive quarks in the spectrum, and we have to understand how their 
appearance can be avoided in an independent particle picture like the HF 
approximation. As already suspected by 't Hooft, the answer lies in the IR 
divergence of the quark self-energy. 



We first remark that the gap equation (135) has no IR problem in the 
limit L —f 00: The Coulomb interaction term has only a simple pole which can 
be regulated with a principal value prescription (or symmetric cut-off). The 
IR divergence discussed in the literature concerns the quark self-energies or, 
in our language, the HF single particle energies which have to be interpreted 
physically as removal energies (Koopmans' theoremlI3). We can deduce them 
by varying the vacuum energy functional 



Ehf/N 



^ ^ (p + 1/2) cos 0{p) - m sin e(p) 
p p 



(139) 



71 7^0, p 
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with respect to the occupation number of level p and find 



^{p) = -j-(P^ -'^Z^) cos6'(p) +msm9{p) + ^ ^ ^ cos(^(p) - 9(j) - n)) 

(140) 

or, in the thermodynamic limit L oo (setting q — 2t:{p + 1/2)/ L), 

• /,/ ^ /■ dg' cos (61(g) - 6'(g')) 
io{q)^q cos 9{q)+m sme{q) + ^ / tT" / a2 • (141) 



277 {q - q'Y 

Due to the double pole, the integral is badly IR divergent. To exhibit the 



divergence, let us isolate the divergent part of the sum in Eq. (140) before 
taking the limit L oo, 

f , ■ nf ,^Ng^ fdq' cos {e{q) - e{q')) - 1 , Ng^L 
UJ(q) — q cos a(q) + msmt)[q) — 



277 {q~q'Y 48 ' 

(142) 

The principal value prescription for the quadratic pole advocated by Bars 
and GreerO and Li et a/.Ej amounts to throwing away the divergent constant 
Ng^L/ 4:8. As a consequence, one gets an awkward sign change (tachyonic 
behavior) in the single particle energies at small p and to. The usual excuse is 
to say that these self-energies have no physical meaning, being gauge dependent 
and non-covariant. Indeed, this difficulty has no direct consequence for the 
vacuum properties — as pointed out above, the gap equation is free of IR 
divergences. Moreover, the calculation of color singlet mesons is also IR safe, 
as pointed out by IL Hooft (his IR cutoff parameter A is related to our box size 
L by A = 12/77L)ej: The infinite self-energy term is cancelled by an equally 
infinite piece in the Coulomb interaction. 

Since the emergence of the constant Ng'^L/ 48 in the single particle energy 



([143) but not in the vacuum energy (139) is a subtle but rather important 



point, let us briefly recall the argument given in Ref. 45 (see also Ref. 74). In 
many-body language, the Coulomb interaction produces both one- and two- 
body terms in H . If the single particle energies are decomposed accordingly 
as 

coip)=J'\p)+Lo(^\p), (143) 
then the HF ground state energy is 

m\0)=-Nj2U'^ (P) + ^^^'^ ip)) ■ (144) 
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As is well known, the factor 1/2 is necessary to avoid double counting of the 
2-body interaction term. In the 't Hooft model, one finds 

27r / 1\ Ng'^L 



u)^^\p) = _ ^ 1^ cos0(p) + msin0(p) 



48 



^^'^W = T^E;Jj(l + ^°^(^(?')-^(?' + "))) ■ (145) 

Remembering that X^n^^^o V"-'^ — '"'^/S, we see that in the ground state energy 
(144) the IR tamed combination (cos(0(p)— 0(p+n)) — l)/n^ appears whereas in 



the quark self-energy ( |143D , the 1-term in w^^^ is cancelled and an IR divergence 
(quadratic pole for L — > oo) survives. 

Why bother about this IR divergence at all? As one proceeds to other 
applications of the 't Hooft model than the meson spectrum, it is no longer 
true that the quark removal energies do not show up in any physical quantity. 
The prime example is finite temperature field theory in th£.large N expansion 
by means of the standard generalization of HF to finite TLJ Here, th£ thermal 
occupation factors depend crucially on the single particle energiesEj If these 
are made finite by some ad-hoc prescription, quarks will manifest themselves 
through contributions to the pressure of order N , in conflict with confinement. 



We shall return to this issue in more detail in Sec. 4.3. Summarizing, we 
believe that the divergence of the quark single particle energies simply refiects 
confinement. In fact this may be the only way how confinement can be realized 
in the independent particle picture characteristic for the leading order in a 1 /N 
expansion. 

3.2.2 Baryons 

Since their structure is very different, we distinguish between those baryons 
which become massless in the chiral limit and those which stay massive. The 
first ones are closely related to Goldstone bosons and will be deferred to Sec. 



3.3| , where they will be discussed for the chiral GN model as well. Here, we 
briefly sketch what is known about the "normal" type of baryons in the 't Hooft 
model. 

Several papers deal with QCD2 baryons fot-a small nurober of colors, like 
N ~ 2 (where they are degenerate with mesonalj) or N = SLlI We concentrate 
on the large N limit instead where baryons are described in relativistic HF 
approximation. An analytic solution like in the GN model seems to be out of 
the question. The work closest in spirit to our description of the GN model 



baryons (cf. Sec. 3.1.2) is Ref. 44. In this work, a numerical approach 
based on a lattice discretization of space and the use of staggered fermions has 
been employed. The HF equation is solved in position space, both at strong 



35 



{Ng^ ^ m) and weak [Ng^ <C m) coupling. One valence level is taken into 
account in addition to the Dirac sea. Unlike in the GN model where the filling 
fraction n/N can be varied, here one has to fill the valence level completely 
with quarks in order to form a color singlet baryon. For heavy quarks, the 
HF equation reduces to the non-relativistic one with a linearly rising potential 
and the sea becomes irrelevant; the HF equation for the valence level then 
becomes 

^ ^/ " yMy)?) vi^) = M^) ■ (146) 

For moderately heavy quarks, significant relativistic effects were found. The 
B = 2 sector was also explored, with the result that two B — 1 baryons were 
observed interacting via an exponentially decreasing repulsive potential. 

3.2.3 Mesons 

Originally, 't Hooft determined the meson spectrum by solving the Bethe- 
Salpeter equation in light-cone gauge, summing up all planar diagrams without 
quark loops. He found an infinite tower of mesons with masses which behave 
asymptotically like 

fi^{n) TT^n + 0{\nn) . (147) 

This should be contrasted to the GN models with one or two mesons only; the 
difference obviously reflects confineme nt ( Regge behavior). In the chiral limit. 



a massless meson appears (see Sec. 3^). These results have subsequently 



been confirmed in normal coordinates, using th^ axial gauge. The necessary 
equations were firaLdsrived by Bars and GreerH and later solved numerically 
by several groupsHa The equal-time approach turns out to be significantly 
more involved than the light-cone approach. Nevertheless, since we cannot use 
light-cone quantization for all the questions which interest us and wish to stay 
within one common framework, we have to outline this approach here. 

As before, we consider small oscillations around the HF solution in the 
space of Slater determinants. The quantized modes of these oscillations are 
just the RPA modes. The RPA equation is identip.1 to the Bethe-Salpeter 
equation in a diagrammatic large N approximationtS To derive it, we return 
to the equation of motion for the bilinear quark operator Q{p,p'), Eq. (|127|), 



insert the expansion (12S) and focus onto the next-to-leading order in l/ViV, 



idtQip',p) - hHFip')Qip',p)-Qip',p)hHFip) (148) 
+ ^^2 XI \{p{p')Qip' -n,p~n)~ Q{p' ~n,p- n)p{p)} . 

n=iO 
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Expanding Q in the HF spinors, as in Sec. 3.1.3 Eq. ( |110| ) and (111), and 
setting 



{n,K\Q2i{p',pm 



then yields the two coupled equations 



(149) 



±E{K)'P±iK,p) 



{ujq{p - K 
Ng^L 



E 



ujo{p))^±iK,p) 
1 



j2 _ 
P'(/P) 



P' 



{s{p,p',K)^^{K,p') 



-c{p,p',K)^±{K,p')} 



(150) 



with the definitions 

s{p,p',K) = sin 



cip,p',K) = 



cos 



dip') - 0{p) 



f 0{p')-Oip) 



sm 



0(p - K)- 9{p' - K) 



f ejp-K) - 9{p' - K) 
V 2 



(151) 



Eqs. (15C) have the standard RPA form of Eq. (116), with 

Ng^L (1 - 6pp 



App>{K) 



Bpp,{K) 



Spp> {ujo{p ~ K) + ujq{p)) 



[P - P' . 



-c{p,p',K) 



Ng^Ljl-Spp,) 
87r2 [p-p'Y 



s(p,p',K) . 



(152) 



The diagrammatic content of these equations is sketchedJii Fig. 5 (see Sec. 
3.2.1 ). Using a shghtly different language. Bars and GreerH interpret them in 
terms of forward and backward going strings. The covariance of the spectrum 
is a non-trivial issue if one works canonically in normal coordinates, but can be 
checked either numerically or via the 1/A^ expansion of the Poincare algebra. 
Meson wavefunctions (or m±her RPA amplitudes) can be found in the literature 
for a few low-lying statesL3 In the limit of heavy quarks, the RPA reduces to 
the non-relativistic Schrodinger equation with a linear potential. The light- 
cone TDA ('t Hooft equation)|-an the other hand can be recovered by going to 
the infinite momentum frameE^I (backward going bubbles are suppressed, see 
also Sec. 3.2). 
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3.3 Chiral limit and Goldstone modes 



Since QCD2 or the chiral GN model in the massless limit break chiral symmetry 
spontaneously (at least at infinite N) , the appearance of a boson whose mass 
behaves like the square root of the bare quark mass has to be expected. This 
has a lot in common with four dimensional theories and makes these models 
good toy models for QCD. However, the chiral limit is somewhat obscured by a 
certain pathology of two dimensional theories which gives rise to a much greater 
variety of (decoupled) massless states. In particular, non-interacting massless 
mesons and baryons already appear at finite N where the Coleman theorem 
certainly forbids talking about Goldstone bosons. Our main interest here is 
in the large N limit where we will derive both massless mesons and baryons. 
However, before doing so, we also briefly explain the more pathological aspects 
of two dimensional theories leading to massless states at finite N as well. 

The identification of massless excitations in two dimensional field theories 

is one of the tasks where light-cone quantization is particularly elegant. Due to 

2 

the light-cone dispersion relation p+ = , a massless state in two dimensions 
has zero energy at all momenta. Therefore its presence can be associated with 
a local symmetry of the Hamiltonian. In the chiral limit th^—right-handed 
fermions drop out of the light-cone Hamiltonian P"*" = H — PLSO'EJ Since the 
Coulomb interaction admits local chiral transformations, P"*" is invariant under 
local phase transformations of right-handed quarks. Free, massless bosons then 
must appear and it is sufficient to identify the winding number of their field 
with baryon number to infer the existence of free, massless baryons as well. 
The IR trouble of Coleman's theorem is avoided by the fact that the massless 
modes decouple. 

The essence of the argument can of course also be phrased without recourse 
to light-cone quantizationE£l Let us return to the fermion currents discussed in 



Sec. 2.2. Vector current conservation allows us to set 

= e^''d,4> (153) 

with a scalar field </>. If the axial current is non- anomalous, we get in the limit 
m ^ (see above, Eq. (pil)) 



a^J^ = 00 = , (154) 

i.e., the announced massless boson. This purely classical argument can easily 
be promoted to the quantum level within the canonical frameworkEj Define 
the right-handed current as 

JrM = ]^{fy{q,t) + jUq,t)) (155) 
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where Heisenberg operators, Fourier transformed with respect to the space 
coordinate x are used. In the chiral hmit, vector and axial vector current 
conservation imply 

[3R{q,t),H]=qjR{q,t) . (156) 
Besides, if P denotes the momentum operator, we have trivially 

[3B{q,t).P]^q]B{q,t) ■ (157) 

For any energy-momentum eigenstate |i) with eigenvalues [Ei, Pi), we thus get 

HjR{qM) = {E^-q)jR{q,m , 

PjR{q,m) = iP,~q)jRiq,t)\i) . (158) 

Thereby decoupled massless bosons can be added to or removed from any 
stationary state. Similarly, left moving massless bosons can be related to the 
left-handed current. 

Let us now turn to the large N limit and discuss the emergence of massless 
particles there. Again we benefit from the experience in many-body theory. 
The Goldstone theorem has a close correspondence which we may invoke here 
(remember that it also holds in non-relativistic theories): If the HF ground 
state breaks a continuous symmetry spontaneously, the RPA develops a gapless 
mode. The basic argument can be given in two lines, using the language of 
linear response theory (without external force) til Take the equation of motion 
for the one-body density matrix, 

ip^[h{p),p] . (159) 

Expand p around its ground state expectation value, p — p^^^+Sp, and linearize 

i(5p= [%(°)),Jp] + 



(0) 

5p 



(160) 

This equation has the form of the RPA, cf. the above derivation in the equa- 



tions of motion approach in Sec. 3.1.3. If the ground state density matrix p 



breaks a continuous symmetry, we find that the deformation Sp corresponding 
to the symmetry transformation (e.g. a chiral rotation) solves the RPA equa- 
tion with zero energy. The backward going ph-bubbles are crucial for getting 
the zero modes; as is well known, the TDA does not have this property. 

These insights enable us to construct the Goldstone boson type solution of 
the RPA equation explicitly. In the case of QCD2, the RPA amplitudes of the 
"pion" can easily be obtained by projecting the right-handed current onto HF 
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spinors. Since these depend only on the Bogohubov angle, the pion-amDhtudes 
are intimately related to the broken symmetry vacuum. One findsuS 



' — cos(0(p)/2)sin(0(p-i^)/2) 



K)/2) 



(161) 



The chiral GN model can be analyzed in a similar manner. Generalizing the 
derivation of the RPA equation, Sec. 3.1.3, to the chirally invariant model, one 
finds an integral equation with a two-term separable kernel (instead of one-term 
in the model with discrete chiral symmetry) and correspondingly two bound 
states, the threshold bound state at 2m (the "cr") and a massless one (the "tt"). 
The pion has the same structure as in Eq. (161), except that the Bogoliubov 
angles now refer to the GN vacuum. They are just the Bogoliubov angles of 
the free, massive fermions given in Eq. ([l^). In this particular case, we are 
in the rare position of having complete analytical control over a meson (the 
tt) in an arbitrary Lorentz frame. In Fig. 0, we illustrate the disappearance 
of the backward components as one boosts the pion to higher momenta (in 
the rest frame, both are exactly equal). At very high momenta, the cpnstant 
wave function familiar from the light-cone meson equation emergesLl In the 
infinite momentum limit, the 't Hooft and GN model pions become identical. 



although their "wave functions" (161) differ in any other frame due to different 
Bogoliubov angles. This illustrates nicely how the vacuum decouples in the 
infinite momentum frame, one of the celebrated advantages of the light-cone 
approach (triviality of the vacuum). Connected to this is the disappearance 
of the backward amplitudes — the RPA reduces to the TDA in the infinite 
momentum frame, the other characteristic benefit on the light-cone. 

We now turn to the issue of baryons which become massless in the chiral 
limit, a phenomenon which does not happen in 3-1-1 dimensions. The best way 
of thinking about them is in terms of coherent, topological excitations of the 
pion field. Like the massive baryons discussed in Sees. 3.1.2 and 3.2.2, they are 
accessible through the relativistic HF approximation. In order to find these 
solutions in the first place, another point p£ view has turned out to be very 
useful. It is closely related to bosonizationPSi but since only the chiral field is 
bosonized, it is conceptually and technically much simpler. _Lt is also similar in 
spirit to the Skyrme pictureupf baryons in 3+1 dimensionsuJ 

Following Salcedo et aLpfl we start from a variational ansatz for the baryon 
near the chiral limit. The interaction part of the ground state energy functional 



E[p\ = N dxtr (-i-f^dx + m'j°)p{x,x) 
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' dxdy\x - y\ti- {p{x,y)p{y,x)) (162) 



(where dx acts only on the first variable of p{x, x)) is invariant under any local 
chiral transformation. Let us perform a local chiral rotation of the vacuum 
density matrix, 

p'ix, y) = e'x(^)^V''(a; - 2/)e-'^(^)'^' . (163) 

We can use the following general decomposition of the baryon density matrices 
p" and p' 

^Pe+ 7°Po - itVi + 7^P5 (164) 

with corresponding superscripts. In the vacuum, pe — S{x — y), due to transla- 
tional invariance all components depend only on x — y and pi vanishes unless 



TO = 0. If the vacuum breaks chiral symmetry, the ansatz (163) yields a new 



density matrix which can be inserted into Eq. (162); x(x) is then determined 



by minimizing E. In order to understand the emergence of baryon number out 
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of winding number, we need the short distance singularities of p" which are 
the same as in the free theory, 

Pq{x) = -^Ka{m\x\) ^ '^\n{m\x\) , 

jn 1 
P5[x) = — sgn(a;)i^i(m|x|) - . (165) 

Using 

p'e{x,y) = cos{x{x) - x{y))pl{x - y) + 2ism{x{x) - x{y))pUx - y) ^ 



p'-,{x,y) = cos{x{x) - x{y))pl{x - y) + ]^ii\^{x{x) - x{y))p"e{x - y) : 



(166) 



one then finds 



Similarly, 



\mi[p'^{x,y)- pl{x'-y)] = -d.j,x{x) , 

y—>x TT 

\imdMy,x)-pUx-y)] = -^{dM^yf ■ (167) 

p',ix,x)Tip[{x,x)=^-f^e^''^^^^ (168) 

where {qq)v refers to the vacuum. If we invoke a finite spatial box of length 
L for these topological considerations, the baryon number gets identified with 
the winding number of the chiral field, 

B = r dx[p',ix, x) - pim = -ixiL) - x(0)) . (169) 
Jo I" 

This is an integer if the baryon goes over into the vacuum at and L. For 
small quark mass, find 

E[p'] = E[p'] +N J da; |i(a,,/))2 + ^cos{V4^c^{x)) - l] | (170) 

with the rescaled field 

(t>{x)^^^x{x) (171) 

and {qq)v referring to the vacuurn-in the chiral limit. E[p'] is minimized if <f) 
satisfies the sine-Gordon equationO 

dl(j> + Vi^m-^^ sm{V^(j){x)) = . (172) 
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The solution with baryon number one is just the famous sine-Gordon kink, 



(j){x) = — ^ arctan exp [ — I . (173) 
It is localized with spatial extent xb-, 

= (-4^™M^y'^' (174) 

and has mass 

If one linearizes the sine-Gordon equation (in the limit m — > 0)p-pne can 
immediately read off the Gell-Mann, Oakes, Renner relation (G0R)l2a 

= -47rm-^ (176) 

for the Goldstone boson. We have thus derived an effective low energy theory 
near the chiral limit. Since the potential energy does not enter Eq. ( [170[ ) 
but only the kinetic energy and mass term, the resulting low energy effective 
theory is universal for a whole class of such chiral models. 

In the limit to — + 0, both the pion and the baryon become massless; the 
baryon gets completely delocalized. This could actually have been anticipated 
from current conservation: di j\ = holds for any stationary state and in two 
dimensions, j\ is just the baryon density jy (cf. Sec. pT^ ). 

An attempt to illustrate the baryon structure in the chiral limit has been 
made in Fig. ^. For to = 0, x(a:) becomes simply a linear function and the 
scalar and pseudoscalar condensates 

{qq) = {qq)vCOs{2TT{x-Xo)/L) , {qij^q) = - (qq) vSm{2Tr{x-xo)/L) , (177) 

can be regarded as projections of a spiral of radius {qq)v winding around once. 
This picture will be generalized later on and turn out to be useful for charac- 
terizing baryonic matter in two dimensional chiral models. 

In Ref. 44 the Skyrme type picture of the baryon was compared to the 
exact numerical solution for light quarks and good agreement was found also 
for small, finite quark masses. It was pointed out that these results support 
the Skyrme picture of the baryon more than the bag picture in two dimensions. 

Finally, let us come back to the question of validity of the variational 



calculation based on the ansatz ( |163| ). In the massless limit, it is easy to 
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Figure 8: Attempt to visualize the topological structure of the massless baryon in terms of 
a "chiral spiral" winding around once. Condensates in units of {qq)v, x in units of L. 



convince oneself that Eq. ( 163 ) with a hncar function x(x) solves the HF 
equation exactly;;. Thus for instance in the 't Hooft model, the massless HF 
equation readso 

--Al^U§^vf\x) + ^ Jdy\x-yKf,{x,y)^<^;\y)^u;^^^-\x) . (178) 
Upon substituting 

^(")(x) = (e-'--^=/^) /^\x) (179) 



as we are instructed to do by Eq. ( |163D , we discover that Lp'-^' does indeed 
solve the HF equation, the only change being that the single particle energy 
LVn gets replaced by uj„ + tt/L. The same argument goes through in the chiral 
GN model, or in any field theory where the interaction term has a local chiral 
invariance for that matter. This shows that the result becomes exact in the 
chiral limit (to leading order in the 1/A'^ expansion, of course). 

4 Finite temperature 

4-1 Finite temperature versus finite extension 

We briefly discussJJie relation between field theories at finite extension and 
finite temperatures Equivalence between relativistic theories at finite exten- 
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sion and finite temperature is formally almost trivial but gives rise to some of 
the most intriguing consequences of covariance. By rotational invariance in the 
Euclidean, the value of the partition function of a system with finite extension 
L in one space direction and [3 in the time direction is invariant under the 
exchange of these two parameters, 

Z(/3,L) = Z(L,/3) , (180) 

provided bosonic (fermionic) fields satisfy periodic (antiperiodic) boundary 
conditions in both time and compact space coordinate. Sending one of the pa- 
rameters L) to infinity, relativistic covariance connects the thermodynamic 
properties of a canonical ensemble with the properties of the pure state of the 
vacuum corresponding to the same physical system but at finite extension. In 



particular, as a consequence of (180), energy density and pressure are related 

by 

e(/3,L)--p(L,/3) . (181) 

For a system of non-interacting particles this equation can be used to relate 
quantitatively the Stefan-Boltzmann law with the Casiniix effect — phenomena 
which at first glance would not seem to be related at allEll Another application 
of this discrete relativistic symmetry worth mentioning -axe quark propagators 
and the interpretation of lattice data in "funny space" EJ In four dimensional 
QCD, the confinement-deconfinement phase transition and the chiral phase 
transition, when quarks are present, appear as "quantum phase transitions" if 
one swaps space and Euclidean time. They are driven by changes in quantum 
rather than thermal fiuctuatiqna which in turn are induced by changes of a 
parameter of the system (L)E3'E3 We have found it useful to apply this tech- 
nique to the two dimensional models at hand to see how it works in detail in 
the presence of interaction effects. Consistency between finite extension and 
thermal field theory calculations provides a non-trivial test of the approxima- 
tions used and will help us to clarify the role of confinement in the 't Hooft 
model. For the case of two dimensions with one extension taken to be infinite, 
the symmetry is illustrated in Fig. ^ 

To leading order in the large N limit, it is easy to include finite temperature 
and chemical potential since the required generalization of our main tool, the 
HE approximation, is well knownL3 Let us focus on finite temperature here 
since finite density and chemical potential will be the subject of Sees. ^ and 
^. The main difference to the temperature zero case is that for sytems in 
thermal equilibrium, one is no longer dealing with a pure state but with a 
mixed state. The concepts of independent particles and single particle states 
remain valid. Whereas at T = each orbit is either empty or occupied (sharp 
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finite L, T = 




Figure 9: Correspondence between relativistic field theories at finite temperature and finite 
extension. 



Fermi surface), at finite T the occupation numh£r is given by the thermal 
occupation probability (Dirac- Fermi distributiorjy ) . The one-body density 
matrix no longer satisfies the projector property — p characteristic for a 
single Slater determinant but assumes the form 



with occupation numbers 



Using 



t + TV.wt , (182) 



one can write 



iV+ + = 1 , N+-N^ = - tanh y , (184) 

P=^ + (7Vo + 7V5)tanh^ . (185) 

where po,P5 are given in terms of Bogoliubov angles in the same way as at 
zero temperature (in the chiral limit, there could be a pi-term in addition.) In 
contrast to the T = case now the single particle energies e become physically 
relevant and will appear explicitly in the self-consistency equation, simply be- 
cause excited states do play a role in the thermal equilibrium. The alternative 
scenario, namely HF at T = but finite extension, does not require any new 
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tools, at least for the GN model. In the case of the 't Hooft model, the situa- 
tion will turn out to be m ore involved due to the different giuon dynamics in 
a finite box (cf. Sec. 4.3). 



4-2 GN model on a cylinder 

Before embarking on the details of the finite temperature GN model, weJiave to 
come back once more to the issue of SSB in lower dimensional systemso Since 
the thermodynamic limit can be taken only with respect to one dimension now, 
there are further restrictions as compared to those of Sec. According to 
the Mermin- Wagner theorem, even a discrete symmetry should be immediately 
restored at T = 0^ . The restoration of chiral symmetry at firute temperature is 
expected to be driven by the presence of kinks and antikinka^jlEJ for any large 
but finite N. If however the limit A'^ — > c» is taken before the thermodynamic 
limit, these configurations are suppressed and one is left with the mean field 
theory where a second order phase transition at a critical temperature Tc ^ 
becomes possible. Similarly, for the continuous model, for large but finite N 
at T ^ 0, the almost long range order disappears. Only if ^ oo is taken 
first, mean field theory applies and one finds the same critical temperature 
as in the discrete model. The non-analytic structure revealed by the phase 
diagram of the infinite A^ GN model depcnda on the diverging number of 
fermion components-like in reduced modelsEjO 

Barducci et alx2 advocate to include a small bare fermion mass, thereby 
eliminating kink-antikink configurations in the thermodynamic limit. Since 
their results go over rather smoothly into those obtained in mean field theory 
at m = 0, we see nothing wrong in considering the chiral limit directly for 
some questions. 

First consider finite extension which requires only minimal modifications 
of the HF approach for the vacuum discussed in Sec. |3.1.l| . Since the fermion 
momenta are discretized due to the antiperiodic boundary conditions, the gap 
equation (|6|) (self-consistency requirement for the scalar density) has to be 
replaced by 

2 Na 

m==^y. , (186) 



2Ng\ y 

n=0 



with 



A^A-^, kn^y{2n + l). (187) 

47r L 



As the finite extension does not affect any UV properties, we can renormalize 
the theory in the limit L oo. Denoting the physical mass on the infinite line 
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by mo and using Eqs. (pq) and (ISC), we obtain 




This equation is now free of divergences but cannot be solved analytically in 
general. To find the critical length below which chiral symmetry gets restored 
it is sufficient to set m ~ and solve for L, 

l-ln—^Y-^. (189) 
2 mo ^ 2n + 1 ^ ' 

n—O 

The sum can be performed, 

J2^—j = ^iC + \nAN^) + 0{l/Ni) (190) 

n=0 

(C « 0.5772 Euler constant) yielding the critical length 

Lc = — . (191) 
mo 

A bare fermion mass mt can easily be taken into account if one replaces the 



gap equation (1S5) by 

r^/^dk 1 1^1 




(192) 

Note the appearance of the new parameter -^p- in addition to mo. Numerical 
solutions of this gap equation for different quark masses are shown in Fig. 
p!o| , exhibiting a second order phase transition in the chiral limit and a cross- 
over phenomenon for non-zero bare masses. Since it explicitly violates chiral 
symmetry, the bare quark mass plays a role analogous to that of an external 
magnetic field in the ferromagnetic phase transition. Finally, we reemphasize 
that to this order in the 1 /N expansion the calculation is identical for the GN 
models with discrete or continuous chiral symmetry. 

We now switch to the equivalent finite temperature calculation to familiar- 
ize ourselves with the corresponding generalization of the HF approximation. 



As indicated in Sec. 4.1, the gap equation of thermal HF has to be modified 



by T-dependent occupation numbers which "smear out" the Fermi surface, 

„ f dk m /3e(fc) , , / — t: 

m-mb^Ng^ / -— tanh , e(fc) = Vm^ + . (193) 

J ZTT e{k) 2 
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T 



Figure 10: Temperature dependence of dynamical fermion mass for several values of bare 
mass, in the GN model. In units of mo- 



The tanh was derived in Eq. ( |l85| ). Renormalizing at T = then yields the 
following equation for the T-dependence of the physical fermion mass, 



, 1 1 

= TOfc 

mo m 



9 I dfc / tanh ^B^/ni? 



27r 



fc2 



(194) 



The integral in (194) converges aniLwe recover the well known results obtained 
by functional integration methods^'L^ in particular a second order phase tran- 
sition at a critical temperature Tc which coincides with from Eq. ( |191| ). 
Bosonic excitations are \/N suppressed and play no role at this stage. Hence 
just like in the finite extension case, there is no distinction yet between GN 
models with continuous and discrete chiral symmetry. 

It is interesting to compare the two approaches leading to the gap equa- 
tions (192) and (194) respectively at an intermediate stage to understand how 
covariance "works" . The crucial identity which ensures the equivalence of the 
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finite L and finite T calculations is 



A/2 



dk 



-A/2 Vm2 + fc2 



; tanh 



(195) 



for (3 = L. Eq. ( |l95| ) can be verified with the help of Cauchy's theorem. The 
integrand has simple poles in the upper half plane at 



k = iVm2 + ((2n + l)7r//3)2 , (n = O.-.TVa) 



(196) 



with residues -2i/(/3^m2 + {{2n + l)7r//5)2). Thereby, a continuous average 
involving thermal occupation numbers can be converted into a sum charac- 
teristic for a finite interval. This is of course reminiscent of the imaginary 
time approach with Matsubara frequencies although different in detail (in the 
canonical approach) due to the interchange of frequencies and momenta. In the 
Euclidean path integral formalism, we would expect no difference whatsoever. 

So far we have concentrated on the dynamical mass or equivalently the 
chiral condensate. Other bulk thermodynamic observables can easily be cal- 
culated from the thermodynamic potential. Let us therefore consider the two 
functions which are minimized with respect to m in the two approaches. At 
finite extension in the chiral limit, it is the HF vacuum energy density 



N 



2Ng^ 



TO^ -|- k?. 



(197) 



(cf. Sec. 3.1.1 ). The gap equation ( 186 ) follows from the condition dE /dm — 0. 
In the finite T case on the other hand, we have to minimize the thermodynamic 
potential ip or free energy density T related to the partition function via 



L 



T ■ 



1 



Intre"''^ 



(198) 



and equal to the negative of the pressure. In the GN model at the HF level, 
the Hamiltonian is given by that of the free, massive theory plus a c-number 
term correcting for double counting of the two-body interaction. 



Lmr 



(199) 



Therefore the partition funiiion can be evaluated following the textbook ap- 
proach for a free Fermi gasp 



T 

N 



2Ng^ 



dk 
2^ 



e(fc) + ^ In ( 1 



-/3e(fc) 



(200) 
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Differentiating with respect to m reproduces Eq. (194) (here for mf, — 0). 
Since we have ah'eady verified that the two gap equations are identical, it is 
sufficient to check that £ and coincide for one value of m (upon identifying 
L and /3) to establish that the two functions agree everywhere as expected 



on general grounds (cf. Eq. (181)). We choose m — where all calculations 
can be done in closed form. For the energy density in the finite interval we 
find, using heat kernel regularization (as usually done in treating the Casimir 
effect), 

£ TT 



N 6L 
For finite temperature Eq. ( |200| ) yields 

T _ 

N ~ ~6P 



(201) 



(202) 



i.e., the Stefan-Boltzmann law. In both cases, the same L- resp. /3-independent 
quadratically divergent term has been dropped. This illustrates the above men- 
tioned correspondence between the Casimir effect and the Stefan-Boltzmann 
law and completes the demonstration that compressed and hot relativistic sys- 
tems are fully related by covariance for the case of the GN model. 

In order to derive other bulk thermodyuamic observables it is useful to 
renormalize the thermodynamic potentialE3l2j Using the gap equation at zero 
temperature one can write 



N 



4, v'°O-0-i^-l/S « 



Following Barducci et alS^ we have added a constant "bag term" mo/47r in 
order to to normalize the (unobservable) vacuum pressure to zero. The deriva- 



tive dJ^rcn/ dm reproduces the renormalized gap equation (194). From this, the 
bulk thermodynamical observables can be obtained. For temperatures below 
the phase transition, pressure and internal energy density are 

e ( (m'^\ \ ml 2 f"" e(fc) 



— = ln^ -l + — + -/ dfc ^-3^ (205) 

N 47r V \ml) J 4n n 1 -|- e^^C^) ^ ' 

where m has to be evaluated at the minimum of J- (solution of the gap equa- 
tion). In the region above Tc, the results are simply 

— = ^ , (206) 

iV 6 47r ' ^ ^ 
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Figure 11: Pressure (solid line) and internal energy density (dashed line) per color in the 
GN model, dividing out the dependence of the Stefan-Boltzmann law. In units of mo; 
adapted from Ref. 92. 



— = \ ^ . (207) 

iV 6 47r ^ ' 

In Fig. |ll|, pressure and internal energy density have been plotted after di- 
viding out the dependence of the Stefan-Boltzmann law, so as to enhance 
interaction effects. 

In conclusion, both finite extension and finite temperature calculations 
agree perfectly for the GN model. Covariance of the single particle energies is 
certainly the key feature behind this result. 



4-3 't Hooft model on a cylinder, decompactification 

The behavior of thaJt Hooft model at finite temperature is rather controversial. 
McLerran and SerO argue that there is naiieconfining phase transition, except 
possibly at infinite temperature. Ming Lifj using standard finite temperature 
field theory methods, concludes that chiral symmetry may get restored in the 
limit T oo. In both of these studies, severe infrared problems were encoun- 
tered, either in the form of divergent diagrams or ambiguous quark self energies. 
Hansson and ZahedEj address the question of the A^-dependence of thermody- 
namic quantities. Several studies have considered QCD2 on a spatial circle at 
zero temperature. As explained above, this can also be reinterpreted as finite 
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temperature calculations for a spatially extended system. Lenz et aim observe 
a chiral phase transition in the massless 't Hooft mod^Lai. some critical length, 
strongly reminiscent of the GN model. Dhar et aLcjtill^d treat the zero mode 
gluons in a more ambitious way than Ref. 45 using technology developed in the 
framework of matrix models and string theory, but are not able to fully solve 
the resulting complicated equations. They propose that the gauge variables 
get decDmpactified by the fermions in complete analogy with the_^chwinger 
modelpj a claim which has recently been disputed by Engclhardtl] 



Let us first consider finite temperature HF. Here, the single particle en- 
ergies enter through the equilibrium thermal occupation numbers. If we used 
the principal value prescription with tachyonic behavior at small momentaEj 



(Sec. 3.2.1 ), we would necessarily get a non-trivial temperature dependence 
in leading order of the 1 /N expansion (T-dependent condensate and pressure 
of order N). Physically, this would imply that the quarks get deconfined. 
Color singlet mesons which might contribute to thermodynamic observables 
can show up only in next-to-leading order in 1/iV. We have argued before that 
the HF single particle energies (interpreted as removal energies) should include 
the constant Ng'^L/48 diverging in the thermodynamic limit (cf. Sec. 3.2.1) 



If this is kept, e ^ and all thermal factors become independent of T. As 
a consequence, 

^liniM = o, lim^^O. (208) 

A system of independent quarks cannot be heated up due to confinement. This 
leaves no room for a chiral phase transition, not even in the limit T — > oo. 

Let us now turn to the finite extension alternative. The 't Hooft model on 
a line reduces to a purely ferniionic theory with linear Coulomb interaction, 



see Sec. 3.2. If we simply take this theory and put it on a finite interval 



with antiperiodic boundary conditions for the quarks, we can repeat the HF 



calculation of the vacuum (Sec. 3.2.1) numerically, in analogy to what we did 
analytically in the GN model (Sec. [l.sj). This calculation has in fact beeii-done 
long ago, although for reasons totally unrelated to finite temperatureca-llHil A 
second order phase transition to a chirally restored phase was found at a critical 
interval length of ~ VdAj^Ng^. Here, the quark self-energies cannot be 
responsible since they do not enter the gap equation which is perfectly well 



behaved in the IR, see Sec. 3.2.1. So something must have gone wrong with 
covariance — the equivalence of finite extension and finite temperature which 
we have verified in the GN model seems to have been l ost in the gauge theory. 



The solution to this puzzle was found only recentlyO£3 It turns out that as 
soon as one considers finite L, it is not true that the gluons in QCD2 can be 
completely eliminated in favor of a static Coulomb potential between quarks. 
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The large N limit of the zero mode gluons (discussed for pure Yang-Mills in 



Sec. 2.5) needs to be refined. Eventually, the finite interval calculation yields 
the same null-result as the thermal HF approach. In effect, the zero mode 
gluons make it impossible to compress the quark system, in the same way as 
the IR-divergent single particle energies make it impossible to heat the system 
up (to leading order in 1 /N) . We are familiar with the fact that gauge theories 
can hide their microscopic degrees of freedom; it is perhaps less well known that 
they can even hide the structure of space-time, making a cylinder look like an 
infinite plane ("decompactification"E£3). Since this is a novel and unexpected 
phenomenon, it may be worth going briefiy through the analysis. 

It is easy to include dynamical fermions into the pure Yang-Mills theory of 
Sec. |2.5| . We again work canonically on a spatial circle of length L ir|-t|bfi| g^uge 
diAi = 0, {Ai)ij = Sijjj^ diagonal in color. The Hamiltonian readai^cJ'EZl 

H^Hg + Hf + Hc , (209) 



with the gauge field kinetic energy 



the quark kinetic energy 



27r 



+ mJ2(4inMn) + bl{n)a,{n)) , (211 



and the Coulomb interaction 



He = - — V ^ ■' — ■ (212) 



The notation for the quark operators and currents is the same as in Sec. 3.2 



As discussed in Sec. 2.5 due to the curved configuration space of the (pi and 
the SU(iV) Haar measure originally appearing in Hg, this Hamiltonian has to 
be supplemented by the following boundary condition for the wavefunctionals, 

"^{ifi, ifN] fermions) =0 if tpi — tpj mod 27r . (213) 

In Ref. 45 quantization was performed after complete classical gauge fixing. 
Thereby, the fact that the (pi are curvilinear coordinates was ignored. This 
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Figure 12: Preferred gluon background field configuration, determining quark boundary 
conditions for 't Hooft model on a spatial circle. Left: Neglecting the Jacobian. Right: 
Taking into account the repulsion due to the Vandermonde determinant ( "pearl necklace" ) . 

led to the assumption that aU the ipi are frozen at the value tt in the large 
N limit. In the resulting purely fermionic theory, the only remnant of the 
giuons are antiperiodic boundary conditions for the quarks in the compact 
space direction. In the meantime, this whole approach has been put on a 
more rigorous basis by first quantizingjji the Weyl gauge and then resolving 
the Gauss law quantum mechanicallycJ This made it clear that the (pi are 
parameters on the group manifold with corresponding Jacobian, the SU(iV) 
reduced Haar measure. When solving the theory, it is then possible to restrict 
oneself to the smallest region in field space bounded by zeros of the Jacobian, 
see Ref. 27 where the consequences for SU(2), SU(3) have been explored in the 
case of adjoint fermions. How can this be generalized to the large N limit? A 
definite choice of "fundamental domain" obviously means that the ipi always 
remain ordered, say < (pi < (^2 < ■■• < 'Pn < Stt. If we think of the giuons as 
particles on a circle, they cannot cross each other and become closely packed 
in the limit N —^ oo. Their fluctuations are suppressed by 1/iV, simply due to 
lack of space. The only degree of freedom left, the collective rotation of this 
"pearl necklace", is a U(l) factor which anyway is not present in the SU(A^) 
theory. This suggests that the correct choice for the gluon background field 
as seen by the fermions is not ipi = tt, but rather the continuum limit of the 
lattice points 

^^ = 27^^, {t^l...N) (214) 

(see Fig. |l^). Instead of antiperiodic boundary conditions, the fermions then 
acquire color dependent boundary conditions which interpolate smoothly be- 
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tween the phases and 2tt, 



V'fc(L)=e'2-'=/^Vfc(0) , {k^l...N). (215) 
In the thermodynamic hmit L oo, both of these choices of the gluon field 



configuration, ipi — it ot Eq. (214), will become indistinguishable and yield 
the well-known results. We now show that at finite L the effects of the gluons 
on the quarks is radically different in these two cases, and that it is the spread 
out distribution ( ^141 ) which is in fact the correct one. 

In such a gluonic background field, the fermions can be treated in a rela- 



tivistic HF approximation along the lines explained in Sec. 3.2.1 except that 



the Bogoliubov angle acquires a color index. The gap equation then becomes 

27r , X . ^ , X ^ , X Q^L sill (0i(n) — 9 j(n ~~ n')) „ . 

— n+a, sine, n -rncos^, n HTTr^y] 7^ ' 216 

L 167r^ (n— aj + ai)'' 

n',j ■> 

where we have switched to the slightly more convenient variable oti = ^ £ [0, 1] 
for the gluons. If = 1/2 as chosen in Ref. 45, 9i{n) becomes i-independent 
and we recover the old gap equation considered in that work. Now, we assume 
ai ~ i/N and perform the large N limit before solving the gap equation. 
Since ai becomes a continuous variable, we replace 9i{n) — > 9a{n) and 

N da', with the result 
2-K 

— (n + a) sin 9a iji) — m cos 9a {n) 
ij 

Ng^L ^ ^^, smi9a{n)-9a,in^n')) ^ ^ ^ ^^17) 



167r2 Jr. (n' ~ a' + a)^ 

n' ^ 

This infinite set of coupled integral equations collapses into a single, one- 
dimensional integral equation, if we set 

9ain) = 9{n + a) . (218) 

Since n is integer and a G [0, 1], this step in effect decompactifies the original 
spatial circle. With this ansatz, the notation n + a = v, n — 11' + a' = 
v' (where v' are dimensionless, continuous variables) and the substitution 
E„' \l da' ^ di/', we obtain 

27r , , , ^ Ng'^L f , sm( 9( ~ 9(iy')) , , 

— lysine - mcos^O^ + Idv' — V a9 = ■ 219 
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The simple pole in the integral has been regularized by the principal value 
prescription. After rescaling the variables via 

27r 27r , , , , 

— -^P , (220) 

where p, p' have the dimension of momenta, we recover exactly the continuum 
version of the HF equation, namely 

P-^-O^^j - moose — 



27r {p - p'Y 

(221) 

Denoting the Bogoliubov angle of the continuum 't Hooft model by 6'cont(p), 
we conclude that 

%) = ^cont f ^i') , (222) 



^ L 

or, in terms of the original, color dependent Bogoliubov angle, 

Q^{n) « flcont [n: + ^) ) ' (A^ ^ oo) . (223) 

This last relation becomes exact in the limit N oo only. In this limit, the 
color- and L-dependent Bogoliubov angles for the 't Hooft model on the circle 
of length L are all given by one universal function, namely the momentum 
dependent Bogoliubov angle of the 't Hooft model on the infinite line. We 
emphasize that this universality only holds for the "pearl necklace" type dis- 



tribution of gauge variables, Eq. (214). If the (fi are all set equal to tt, there 



is no analytic way known how to relate 9{n) for different L values, but one has 
to solve the gap equation numerically for each LeB 

The upshot of this simple exercise is the following: In the large N limit, 
the gluon variables influence the fermion boundary conditions in such a way 
that the circle gets replaced by a line — they decompactify space-time (see 
Fig. this only works for gauge invariant variables, where one sums over all 
colors). The length L of the spatial circle becomes an irrelevant parameter. 
To confirm this last point, let us evaluate the condensate as an example of an 
observable, 

im) = -—^sm9i{n) = --^X! J dQ;sin0(n-t-a) = -N J — sin 6lcont (p) ■ 



(224) 

The sum over the discrete momenta and the color sum in the large N limit con- 
spire to produce the continuum result, independently of the starting L value. 
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Figure 13: Sketch (in momentum space) of how color-dependent quark boundary condi- 
tions in a finite interval may simulate a theory on the infinite line, in the limit Af — > oo 
( "decompactification" ) . 



In the alternative thermodynamic interpretation, the condensate becomes in- 
dependent of temperature to leading order in 1/N, now in perfect agreement 
with the finite temperature HF result. 

Notice that so far, we have discussed the influence of the gauge fields on the 
quarks which is indeed dramatic. Vice versa, we do not expect the quarks to 
influence signiflcantly the gluon zero point motion, again due to the effects of 
the Jacobian. The kinetic energy Hg will give the same result as in pure Yang 
Mills theory. Since this contribution to the energy density is L-independent, it 
again yields zero pressure, reflecting the absence of physical gluonic excitations 
in 1+1 dimensions. 

4-4 1/-^ corrections at finite T, role of pions 

In the present section, we wish to discuss the role of the pions at low tem- 
perature. Since this has to do with the next order in the 1/N expansion and 
cannot be covered systematically here, we will content ourselves with briefly 
reviewing the state of the art and indicating the issues which have attracted 
most attention recently. 
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The importance of pions_fat_low temperature QCD was pointed out by 



Leutwyler and coUaboratoralHjilHj some time ago. In the vicinity of the chiral 



Umit, pions are naturally the lightest bosons and as such preferred degrees 
of freedom at low temperature. Using chiral perturbation theory, interaction 
effects in the pion gas can be computed reliably to rather high order. An ex- 
trapolation of the temperature dependence of the chiral condensate even gives 
a surprisingly good prediction for the critical temperature of the QCD chiral 
phase transition. Glueballs (if they exist in nature) are much too heavy to 
matter at low temperatures, quarks are anyway ruled out due to confinement. 

The situation is somewhat different in the case of the NJL model which 
does not exhibit confinement. However also here thermal excitation of a qq pair 
is strongly suppressed as compared to pions at low T. It was soon recognized 
that the popular expansion is somewhat misleading physically: To leading 
order, the pions do not contribute at all to thermodynamic observables, they 
are dominated by the quarks due to their color weight factor (entropy). For any 
finite N, one can estimate a temperature below which pions dominate, hut.|for| 
N ^ oo this temperature moves to 0. Zhuang, Klevansky and Hiifner^°^ 
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and others have developed techniques to go to next-to-leading order in the 
1/N expansion, e.g. in the thermodynamic potential. They take into account 
fluctuations about the mean field which correspond diagrammatically to a 
ring sum and physically to mesonic contributions. Both qq bound states and 
scattering phase shifts enter. Since a meson in the NJL model can decay into 
quarks, a "mixed phase" arises at intermediate temperatures where mesons as 
wqlJ_^iS_fluarks play a role. This problem was also investigated by Barducci et 
including the temperature dependence of meson masses. 
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The same issues are relevant for the lower dimensional models which are 
the subject of the present work as well. Since these are toy models anyway, one 
usually takes the limit N oo seriously (not as an approximation of = 3), 
at least this is the point of view which we have adopted here. In this case, one 
might argue that the pions have no influence on the phase diagram of the chiral 
GN model. In the strict chiral limit, this is a rather extreme picture: N has 
to be so large that massive fermions win over massless pions in spite of their 
exponential suppression through the Boltzmann factors, but as a mathematical 
limit it is perhaps acceptable. 

It is instructive to think about the corrections (meson effects) in the 
"rotated" picture of finite extension field theory. There, in next-to-leading or- 
der one would have to evaluate the RPA vacuum in 



finite interval, in clear 

correspondence to the bubble sum of Zhuang et alS£B Unlike the thermody- 
namic calculation, this does not require any new formalism since the RPA 
technique is well established. Fermionic and mesonic degrees of freedom are 
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then treated consistently and it is also clear that scattering states as well as 
bound states enter. Such a calculation has not yet been carried out for the 
GN model. The thermodynamic treatment of the corrections has been 
presented by Barducci et aZ.E2j showing low temperature dominance of pion- 
like excitations in the massive GN model at finite N. Introducing a bare 
quark mass enables the authors to vary the degree of explicit chiral symmetry 
breaking. 

Let us finally come back to the 't Hooft model. Here, due to confinement, 
the quarks do not contribute at all to the pressure at order N so that the 
mesonic contributions are in fact the leading ones. Since in the large N limit 
one has non- interacting mesons, one can use the thermodynamics of the ideal 
Bose gas to analyze this situation. Near the chiral limit, the pion will naturally 
stick out at low temperatures. 

A theoretical question of some interest is the corresponding finite inter- 
val RPA calculation in the 't Hooft model. Here, we do not expect that the 
mechanism which leads to "decompactification" in leading order remains in- 
tact, otherwise one would also not get any meson contribution to the pressure. 
This must mean that the fluctuations of the zero mode gluons come into the 
picture. It would be interesting to study this in more detail since it would teach 
us something about the gluon dynamics at large but finit&JV and perhaps help 
to settle some of the difficult issues raised by Dhar et h/.EjOcS in connection 
with the string picture of QCD2. 



5 Finite density 

In the present section, we address the problem of finite baryon density or 
chemical potential in the GN and 't Hooft models. Since the physics issues are 
quite different from those at finite temperature (Sec. ^), we restrict ourselves 
to the T = case here and try to clarify the structure of cold baryonic matter 
first. Ultimately, we are of course interested in the phase diagram of our 
soluble models in the whole (T, /i)-plane; this will be the subject of Sec. 6. 
Whereas the GN mpHp] at. finite density has received a lot of attention in 
the literaturepZl'Ooliij'LJ the corresponding extension p£,ihe 't Hooft model 
has remained essentially unexplored until very recentlyoLJ Throughout this 
section, our main focus will be on the relationship between the structure of the 
single baryon and that of baryonic matter. 



5. 1 GN model as a Fermi gas 

To determine the ground state of the GN model at finite density, we can again 
use the HF method expected to become exact in the large N limit. As we shall 
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see, this does not mean that one should trust the results of a HF calculation 
blindly: If the ground state breaks some symmetry, the HF equations will 
develop different solutions and one may easily miss the solution with the lowest 
energy due to erroneous assumptions. Anyway, to start our investigation let 
us first assume translational invariance and only offer the system the chance 
to break chiral symmetry. We thus simply fill a number of positive energy 
plane wave states on top of the Dirac sea until we reach the desired baryon 
density. This is the standard approach to the finite density GN model, phrased 
in many-body language. Although we have in mind primarily the chiral GN 
model, the corresponding calculation would in fact be identical for the model 
with discrete symmetry only. 

We denote the fermion density per color (or baryon density) by ps — pf /it 
{pf: Fermi momentum). At the mean field level, the fermions acquire a physical 
mass m which has to be determined self-consistently. The ground state energy 
density per color is then given by 



N 



2-K 



Pf 



2Ng^ 



(225) 



where A is the UV cutoff. The first term sums up the single particle energies for 
all occupied states (the Dirac sea plus all positive energy states with \p\ < Pf), 
the second term corrects for double counting of the interaction. We renormalize 
the theory at pf = 0, cf. Sec. 3.1.1, and denote the physical fermion mass in 
the vacuum by mg. Using the vacuum gap equation ( |69| ) to renormalize the 
matter ground state energy density, Eq. ( |225D , we find (dropping an irrelevant 
infinity) 



£_ 

N 




(226) 



The energy is minimal provided m satisfies 



m In 



m? + p^j 



mo 



(227) 



i.e., for 



m — or 



mo 



(228) 



61 



Pi 

Figure 14: Physical fermion mass as a func- 
tion of Fermi momentum in the GN model. 
In units of mo. 



Ps 

Figure 15: Energy density per color as a 
function of Fermi momentum in the GN 
model. Solid curve: Chirally symmetric 
solution (m = 0). Dotted curve: Broken 
chiral symmetry (m according to Fig. 14). 
Dash-dotted line: Mixed phase. In units of 
mo. 



The corresponding energy densities are 
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(229) 



The physical quark masses (228) and the energy densities ( |229|) are plotted 
in Figs. |l^ and From these figures one might be tempted to conclude that 
chiral symmetry is broken at low densities and gets restored in a second order 
phase transition at pf — mo/2. As is well known, this does not occur, rather 
there is a first order chiral phase transition a,t pj — mo/\/2- This can easibt^be 
inferred by inspection of the thermodynamic potential of the GN modelEjO 
At zero temperature, it is possible to understand it also more dir ectly in terms 
of the following instability. Let us compare the energy densities ( |229[) with the 
energy density for a system of size L divided into two homogeneous regions 
I (size i) and II (size L — £). In region I chiral symmetry is restored; it 
contains the-epctra fermions needed to get the prescribed average density (the 
"MIT bag"lliil). Region II consists of the physical vacuum with broken chiral 
symmetry, void of excess fermions. The mean energy density obtained in this 
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way IS 



(230) 



(231) 



N \ L ) Att 2-kI 
Minimization with respect to i yields 

Too 

valid for p/ < toq / \/2, and hence the optimal energy density 

iV 47r V27r V ^2/ ^ 

As shown in Fig. |l^, this solution is lower in energy than the homogeneous one; 
moreover, it yields the convex hull off. It ends exactly at the first order phase 
transition point p/ = too/V2 where all space is filled with one big bag. This 
should be contrasted to the scenario underlying Fig. |lj where the fermion 
mass decreases continuously. We thus recover the well known mixed phase 
interpretation of the GN model at finite density, in the zero temperature limit. 
Notice also that only the total size of regions I and II nmtters, not how they 
are subdivided; there could be baryon "droplets" as wellti3 

One important point to which we would like to draw the attention of the 
reader is the behavior of £ near ps — pf /tt — 0. Since ultimately, at very low 
density, the fermionic matter problem must reduce to the problem of a single 
baryon, one would expect 



d£ 



dpi 



= Mb (233) 



where Mb is the baryon mass. In the present calculation, AIb is not the 
physical baryon mass, but the mass of an alleged "delocalized" baryon. This 
is inherent in the translationally invariant HF approach, i.e., the assumption 



that the single particle orbitals are momentum eigenstates. Using Eq. (233) we 



obtain in the homogeneous, single phase calculation, Eq. (229), Mb = Nttiq 



consistent with a short range force and a delocalized baryon. The (physically 



more viable) mixed phase approach, Eq. (232), predicts a baryon mass lower 
by a factoE_Df 1 / \/2 which can be understood in terms of the bag model for 
the baryonU 

However, the GN model possesses bound baryons with lowest mass Nrnojji 



(kink solution for the model with discrete chiral symmetry, cf. Sec. 3.1.2) 



or even massless baryons (model with continuous chiral symmetry, cf. Sec. 
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3.3). These binding effects are not suppressed and should be correctly 
reproduced in a HF approach, in the low density limit. They have obviously 
been missed here due to our tacit assumption of translational invariance. There 
is no good reason why such effects should not play a role at higher densities 
as well. Moreover, differences between the continuous and discrete chirally 
symmetric GN models based on their different baryon structure and spectra 
are not at all captured by the Fermi gas approach. In Sec. we shall 
present a cure for this disease. Before that however, let us first repeat the 
naive calculation for the 't Hooft model where we can also investigate the role 
of confinement. 

5.2 't Hooft model as a Fermi gas 



We have discussed the vacuum and baryons in the 't Hooft model in Sec. 3.2 
Assuming translational invariance, it is again straightforward to include a finite 
baryon density into the HF calculation. If p f d enotes the Fermi momentum 



we only need to replace the density matrix (134) by 



P{p) = ^{Pf - \p\)u{pW{p) + v(jp)v'^{p) 

= Q(j>f-\P\) + Q{\p\-PMPW{P) (234) 

where we have used the completeness relation for the spinors in the second 
step. In the expression for the HF ground state energy density at m = (|139|), 



according to the second line of Eq. (234), we must exclude the region [—pf,pf] 



from the momentum integrations and pick up an additional term due to the 
change in the baryon density trp. 



N 



-;^Q{\p\~ Pf)p cos e{p) 



This yields at once the following finite density generalization of the gap equa- 
tion, 

psin.(,) + ^ - (W>^/)' (236) 

whereas the condensate now becomes 

{qq) ^-N I '^e{\p\ - Pf) Sineip) . (237) 
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Figure 16: Quark condensate as a function of Fermi momentum in the 't Hooft model, in 
units where Ng^ = 2n. 



The gap equation ( 236 ) can easily be solved numerically for various p / . The re- 
sulting condensate is shown in Fig. ^ We find that it decreases monotonically 
with increasing density, disappearing at a critical Fermi momentum 

P/- 0.117 (^y^' ■ ^238) 

This behavior is strikingly similar to the corresponding result for the GN model 
depicted in Fig. |lj, again suggesting some phase transition with restoration of 
chiral symmetry at high density. Unfortunately, here we are not able to assess 
whether we are dealing with a first or second order phase transition. The 



reason lies in the following problem: If we compute the energy density (235) 
for the 't Hooft model we discover that subtraction of the value at p/ = is 
not sufficient to give a finite result. Unlike in the GN model, the difference is 
still IR divergent. To be able to proceed, we enclose the system in a box of 



length L. We then find that the divergence is due to the last term in Eq. (235) 
(the one which does not involve the Bogoliubov angles) which now contributes 
the following double sum to the energy per color. 



E 
N 



Ng^L 



E E ^- (239) 
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Here antiperiodic boundary conditions for fcrmions have been employed in 
the box regularization, and correspondingly the interval / is defined in the 
following way, 

/=[— n/,n/] for B = 2ri/ + l odd , I~[—nf — l,nf] ioi B — 2n j +2 even . 

(240) 



The result (239) is even more alarming than the non-convex behavior of £ 
in the GN model, Fig. ^ due to its L-dependence. Adding quarks to the 
vacuum causes the energy to increase by an infinite amount in the limit L 



oo. Evaluating the double sums in Eq. (239) for low values of B, we obtain 
information on the origin of this divergent behavior. For B = 1 (/ = [0, 0]) in 
particular, the calculated baryon mass (to leading order in L) is 

Mb = N f . (241) 



V 48 

This is the same relation as AIb = NnriQ in the GN model except that the phys- 
ical fermion mass is replaced by the infinite constant Ng^L/AS characteristic 



of confinement, cf. Eq. (142) 



Summarizing, the problems encountered in the GN model with translation- 
ally invariant baryonic matter again show up in the 't Hooft model, although 
in aggravated form. The physics reason is clear: In the GN model the cost 
of distributing TV fermions over the whole space is governed by their physical 
mass; in the 't Hooft model, due to confinement of quarks, the correspond- 
ing quark effective mass diverges with the volume. On the other hand, it is 
known that both models do possess massless, delocalized baryons in the chiral 
limit. Evidently, this has to be accounted for, and we conclude that the naive, 
translationally invariant HF approximation cannot be trusted. 

5.3 Breakdown of translational invariance 
5.3.1 Strict chiral limit 

We now return to the massless baryons introduced in Sec. 3.3 in the chiral 
limit of two-dimensional field theories. By using literally the same techniques 
(except for a different value of the baryon number) we can easily find the 



ground state of the system for any baryon density. As discussed in Sec. 3.3 
and illustrated in Fig. ^, the single baryon can be visualized in terms of a 
"chiral spiral" parametrized by x(^) with one single turn in the whole volume. 
A finite density ps = B / L = pf /t: on the other hand implies that 

X{x) = pf{x - xo) , (242) 
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Figure 17: Same as Fig. |8| but for baryonic matter in the chiral GN or 't Hooft model 
("chiral spiral"). The wavenumber of the oscillations in the condensates is 2pf. 



i.e., one full rotation over a physical distance which now has a well defined 
limit for L oo^ namely 2/ pg. The baryon density remains constant in space 
for the reasons discussed in Sec. 3.3, but the condensates are modulated as 



{m) ^ {qq)v cos 2p f {x - xq) , 

(91759) = - {qq)v sin 2p f {x ~ xo) 



(243) 



They can again be viewed as projections of a chiral spiral of radius | {qq)v \ onto 
two orthogonal planes, see Fig. 17. 



This state breaks translational symmetry; it is a crystal. In fact, it may be 
viewed as the simplest possible realization of the old idea of a Skyrme crystal^ij 
here in the context of large N two-dimensional field theories. It is meaningless 
to ask where one baryon begins and ends; each full turn of the spiral adds 
one unit to the baryon number. Only the condensates reveal that translational 
symmetry has been broken down to a discrete subgroup. The energy density of 
this unusual kind of "nuclear matter" is simply (after subtracting the vacuum 
energy density) 

^ ' ' (244) 



__p1 

N 2tt 



Surprisingly, this is exactly what one would get for a free Fermi gas of massless 
quarks although Eq. (244) holds for interacting theories where the vacuum has 
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Figure 18: Same as Fig. ^ (GN model), except that the energy density of the Skyrme 
crystal type of state (the true ground state) has been included as the dashed line. 

lower energy due to chiral symmetry breaking. In Fig. ^ we compare the en- 
ergy density for this state to the ones discussed above for the GN model, where 
translational symmetry had been assumed. The crystal is always energetically 
favored, the dependence on pf is now convex, and there is no trace of a phase 
transition, neither first nor second order, at any density. The horizontal slope 
a.t Pf — correctly signals the presence of massless baryons and eliminates 
the above-mentioned problems with the spurious massive, delocalized baryons. 
We cannot even draw the corresponding picture for the 't Hooft model, simply 
because in this case the quark Fermi gas is infinitely higher in energy than the 
Skyrme crystal for L ^ oo. Nevertheless, all the results for baryonic matter 
discussed in this section apply to the 't Hooft model as well. 

In the high density limit the oscillations of the condensates become more 
and more rapid. If we are interested only in length scales large as compared 
to 1/pf, the condensates average to zero. In this sense, one might argue that 
chiral symmetry gets restored at high density, although not in the naive way 
suggested by Figs. |lj or |l5|. 

Finally, we remark that the "chiral spiral" ground state for fixed baryon 
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density still preserves one continuous, unbroken symmetry, namely the combi- 
nation of translation and chiral rotation generated by P+PfQs {P: momentum 
operator, Q^: axial charge). One would therefore expect that RPA excitations 
on this ground stateCjS (or mesons in nuclear matter) will have only one 
collective, gapless mode, a hybrid of a "phonon" and a "pion" . 
5.3.2 Non-vanishing bare quark masses 



As discussed in Sec. 3.3 following Ref. 44, a finite bare quark mass can be 
included in the Skyrme type approach to the baryon; it changes the effective 
action for the chiral phase x{x) from that of a free massless field into a sine- 
Gordon model. The baryon was identified with the kink solution. It is then 
clear which type of solution is likely to be a good candidate for baryonic matter: 
The sine-Gordon kink crystal. It goes over into the results of the previous 
section in the limit m — > and is an approximate solution to the HF equations 
(presumably a very accurate one, judging from the results of Ref. 44 for the 
single baryon). Luckily, the sine-Gordon kink crystal has, already been studied 
thoroughly in the literature, first in solid^tate physicstliHlia and more recently 
as a toy model for the-Skyrme crystaipia in terms of Jacobi elliptic functions 
and elliptic integralstiJ We take over the results from Ref. 116 which is close 
in spirit to the present study although the authors apparently did not have in 
mind two-dimensional large N field theories. Adapting the formulae of their 
work to our notation, the following steps allow us to generalize the Skyrme 
crystal of the previous section to small, finite bare quark masses: Let 
denote the mass of the Goldstone boson, Eq. ( |176| ), and ps = Pf/ir the average 
baryon density (this is our definition of for the case of broken translational 
symmetry). We then first have to solve the transcendental equation 

2fcK(fc) (245) 



Pf 



for k where K(fc) is the complete elliptic integral of the first kind. The sine- 
Gordon kink crystal is then given by the following solution of Eq. (|172[) , 



X(a;)--+am(e,fc) , ^ ^ ^{x - xo) , (246) 

(am(^. A;) is the Jacobian elliptic amplitude function). From this, we can ex- 
press the baryon density and the various condensates in terms of further Ja- 
cobian elliptic functions (dn, sn, cn) as follows, 

psix) = -dxXix) = ^dn(^,fc), 

(qq) = {qq)vCOs2x{x) = ~{qq)y {cii'^{^,k) ~ sii^{^,k)) , 
(91759) = -(qq)v sin 2xix) ^ +{qq)y2sniC,k)cn{^,k) . (247) 
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Figure 19: Solid curve: Spatial oscillation of the baryon density in the regime Pf JTItt 
Circles: Baryon density for a single baryon for comparison. 



Here ^ is as defined in Eq. (246). Finally, the energy per volume of this kind 
of matter is given by 

E(A;) denoting the complete elliptic integral of the second kind. 

Let us now illustrate these results in two regimes of interest, namely at 



low and high density. At low density (p/ <C mTr), k in Eq. (245) approaches 1 
exponentially, and the baryon density features a chain of well resolved lumps 
whose shape is determined by the single kink solution (Fig. Likewise, 
the condensates behave like those of a single baryon: {qq) changes from the 
vacuum value outside the baryons to its negative in their center whereas {q I'y^q) 
is peaked in the surface region of each baryon (Figs. |2^). These condensates 
are projections of the distorted "chiral spiral" shown in Fig. G2. The energy 



(248) for low densities behaves as 



E ^ = MbPb , (249) 

showing the expected connection to the baryon mass. At high densities {pf ^ 
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Figure 20: Solid curve: Spatial oscillation 
of the scalar chiral condensate in the regime 
Pf ^ mir- Circles: Scalar chiral condensate 
for a single baryon. 



Figure 21: Same as Fig. 20, but for the 
pseudo-scalar chiral condensate. 



m^), k approaches like 

fc « ^ . (250) 
PI 



Thus ^ in Eq. ( |246| ) becomes pf{x — xo)- Moreover, for k ^ 0, the Jacobian el- 
hptic functions am(^, fc), sn(^, k), cn(^, k) are known to reduce to the argument 
^ and the ordinary trigonometric functions sin^ and cos^, respectively. We 



thus recover the results for the simple chiral spiral in Sec. 3.3 (the parameter 
xq has to be readjusted to take care of the shift by tt/2 in Eq. ( |246| )). The 
energy in this case is approximately 

TV ^ 2^ + Stt ■ ^^^^> 

The condensates look very much like the sin- and cos-functions of the massless 
case and need not be plotted. The baryon density oscillates around a constant 
value, reflecting the strong overlap of the baryons, and are well approximated 
at high density by 

PB(:r:)«^(^l-^(^)'sin2p/(x-xo)^ . (252) 

The behavior of the baryon density pb{x) as one increases pf (i.e., the mean 
density) is illustrated in Fig. In the chiral or high-density limit {171-^ /pf —> 
0) pb{x) eventually becomes a;- independent. This provides us with another 
way of understanding the structure of matter described in the previous section, 
namely as arising from a chain of very extended, strongly overlapping lumps. 
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Figure 22: Illustration of the distorted "chiral spiral" for baryonic matter at non-zero bare 
quark mass (sine-Gordon kink crystal). 



5.3.3 Comparison to other works 

It is noteworthy that a similar chiral structure of fcrmionic matter has been 
reported previously in a variety of models different from the ones discussed 
here. This indicates that the basic results are more generally valid than our 
derivation might suggest. Let us briefly go through these works to see what is 
common and what is different from our case. 

The first mention of nonuniform structures we are aware of is Ref. 118 on 
the massive Schwinger model. Since in the standard Schwinger model baryon 
charge is confined, one either has to add an inert, uniform background charge 
to neutralize the system or introduce flavor. In the first case, Fischler et 
al. find that a spatially varying, periodic charge Jiinsity is induced which 



they interpret as 1-1-1 dimensional Wigner crystalti3£ify More recently, the 
massless Schwinger model at finite chemical potential has also been examined 
and nscillati pg chiral condensates with wavenurnhfir k = 2fi were found, by 
bosonizationE£iJ as well as by functional methodsE£3 In the massless case, the 
charge density becomes uniform. This phenomenon has obviously a lot in 
common with our findings, although it is different in detail. In particular, 
the way to avoid conflict with the no-go theorem for spontaneous symmetry 
breaking in two dimensions must be different since one cannot invoke large 
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Figure 23: Spatial dependence of the baryon density as it evolves with increasing average 
density (or Fermi momentum), in units of m,r- 



iV arguments; presumably, in the U(l) case, the anomaly and the long range 
Coulomb interaction play a crucial role. In the generalized Schwinger model 
with two flavacs, one can have baryonic charge for neutral systems. Here, 



Fischler et alxlB also arrive at a sine-Gordon equation for the light meson but 
point out that a possible crystal solution would be destroyed immediately by 
quantum fluctuations. This is of course avoided in our case by the large iV 
limit. 

QCD2 witi. flavor at finite chemical potential has been discussed by Chris- 
tiansen et al&23 in the chiral limit. Oscillatory condensates where found within 
the path-integral approach. Since these authors work at finite iV, it is not clear 
how the no-go theorem is avoided. On the other hand, since we know that 
QCD2 does have massless mesons and baryons at finite N, one cannot rule out 
that there is yet another way of bypassing the Coleman theorem. 

Even more surprising are perhaps quite a number of reports about spa- 
tially inhomogeneous chiral condensates with exactly thp same wave number 



as in our case, but in 3+1 dimensions. Kutschera et allEz3 have studied quark 
matter with pion condensates within the cr-model, using a mean-field approx- 
imation. Whereas their quark matter is closely analogous to our Fermi gas, 
the pion condensed state resembles the chiral spiral, especially as far as the 
condeaaates are concerned. Another important work is the one by Deryagin 
et aZJlHa on large N QCD in 3+1 dimensions. With the help of a variational 
calculation, these authors find an instability of cold and dense matter with 
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respect to formation of an inhomogcnous, anisotropic condensate, a standing 
wave with wave number 2/i. They point out that unhke color superconductiv- 
ity, this condensation survives in the N oo hmit. PhysicaUy, they interpret 
their resuhs as condensation of pairs of particles and holes on opposite points 
of the Fermi surface \p\ = fj, ( "Overhauser effect" E£3). This scenario was sub- 
sequently taken up by Shuster et aLE£D and Park et aZJlHs who elaborated on 
the competition between BCS (pp or hh) pairing and Overhauser (ph) pairing. 
They confirm that the large N limit favors the chiral wave, although it seems 
very unlikely to be relevant for N values as low as 3. 

Let us see how this pairing picture fits into the 1+1 dimensional models 
considered here. Since relativistic HF is very similar to BCS theory, we can 
interpret the vacuum result (as encoded in the Bogoliubov angles, cf. Fig. |[ 
Sec. |3.2.l[ ) as a ph-pairing in a momentum region with a width of order m 
or y/Ng^ (for the GN and 't Hooft models, respectively) and centered around 
p = Q. The crystal solution is obtained by substituting 

q{x) e'Pf'^'qix) , (253) 

thereby splitting momenta between right-handed and left-handed fermion com- 
ponents by ±p/. Then, provided that Pf ^ {m, y'Ng^), we can indeed de- 
scribe the result as ph pairing on opposite sides of the Fermi "sphere" (here 
an interval). However our results also hold for small pf where this picture is 
not really applicable. 

6 Phase diagram in the (T, /i)-plane 

Most of the work done so far on the phase diagram for the two dimensional 
large N models has been devoted to the original, non-chiral GN model. Its 
amazingly rich phase diagram has lieen discussed comprehensively by WolfiE3 
(for rrih = 0) and by Barducci et alx3 (for finite bare quark masses). It has also 
been invoked foiLJinderstanding some questions about the statistical physics 



of polyacetylenel^zil a linear polymer. We are primarily interested in the role 
played by chiral symmetry and will therefore concentrate on the chiral GN 
model. Actually, as far as the conventional Fermi gas approach is concerned, 
both variants of the GN model yield identical results to leading order in l/iVE3 
Apart from briefly reviewing the common lore about the GN model, we will 



also investigate the novel "crystalline" phase (chiral spiral, cf. Sec. 5_^) and 
determine which phase is thermodynamically stable. In order to keep things 
technically as simple as possible, we consider only the mf, = case and work to 
leading order in the 1/N expansion throughout this section. Incidentally, the 
't Hooft model will be no issue here: Since in the limit N oo confinement 
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suppresses any T-dependence (cf. Sec. 4.3), there is nothing we could add to 
what has already been said about T = and finite density in Sec. ^ 

6. 1 Chiral GN model with translational invariance 

Let us assume unbroken translational invariance and, without loss of generality, 
{qij^q) = 0. If we are interested in equilibrium thermodynamics at finite 
temperature and density, it is convenient to work with the grand canonical 
ensemble where the baryon density is adjusted via a Lagrange multiplier, the 
chemical potential /i. The HF method can be generalized without any difficulty 
to finite ^LJ Thereby, the HF equation remains unchanged and equal to the 
free Dirac equation, 

7^T^ + n^7°^ V^n{x) = en(Pn{x) , (254) 

whereas the chemical potential only enters the gap equation through the stan- 
dard Fermion occupation numbers 

rn_{qq)_^__ 1 



Ng2- N -^^"(^V»(^) e/3(.„-,0 + l ■ (255) 

The solution of the HF equation is once more given by the free, massive Dirac 
theory with e„ = ±e(fc) = ^/W^+rn? . The relevant thermodynamic function 
here is the grand canonical potential often referred to as "effective potential" 
in this context, 

l^eff = - ^ In ^ = - ^ In tr e^^^^"^^) . (256) 

Recalling that the Hamiltonian of the GN model differs from that of the free 
massive Dirac theory only by a c- number (in leading order in 1/A^) and using 
the partition function for a free Fermi gas as follows, 

InZ = ln]j(l + c-'^('"-'')) (257) 

11 



2tt 

the effective potential density (Vcfr = Vcs/L) per color becomes 



, , , e(k) + u + - In 



(258) 
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In this way, the "vacuum part" [T = = 0) can easily be identified and 
separated from the "matter part" which in turn is made up of particle and 
anti-particle contributions. Incidentally, the term 



dfc 
2^ 



(259) 



in eq. ( ^58[ ) is a pure vacuum term which is usually omitted in the literature; 
it arises from the infinite fermion density of the Dirac sea. The gap equation 
(255) can be further worked out to yield 



2tt 



dk 



e{k) 
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1 



,/3(e(fc)-^t) 



g/3(e(fe)+p) _^ I 



(260) 



where we have used 



u{k)u(k) = —v(k)v{k) 



(261) 



Eq. (260) is equivalent to the condition that the effective potential ( p58| ) 
stationary with respect to m. 



aVcff 
dm 



= 



(262) 



Renormalization can be performed at T = 0, /i = (physical fermion mass mo) 
resulting in the renormalized potential 
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(263) 



)(> 



-/3(e(fc)+p) 



By adding the "bag pressure" mQ/(47r) and dropping the term in Eq. (259), 
it has been adjusted to vanish in the vacuum. The effective potential ( 263 ) 
contains all the information about the conventional phase diagram of the GN 
model. (For the path integral derivation and the representation of Vcff in 
terms of sums over Matsubara frequencies instead of integrals over thermal 
occupation numbers, we refer to Refs. 95, 110.) For the sake of completeness, 
we record the expressions for standard thermodynamical bulk observables: The 
free energy density J-^ (or, equivalently, minus the pressure p) is equal to VeS, 
Eq. (263). The quark density p and the entropy density s are obtained by 
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Figure 24: Phase diagram of the (chiral or non-chiral) GN pnH els. assuming unbroken 
translational symmetry. For a discussion, see text below Eq. (266). Units of mo; adapted 
from Ref. 95. 



differentiating p with respect to fi and T, 
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g/3(£(fc)-M) -I- 1 e^('^(''')+^) + 1 
e(fc) — /i e(fc) + [I 



71" .In \ e''(^('')"^) + 1 e''('^('')+^' + 1 
whereas the internal energy density can now be deduced via 

e = Ts — p + jip . 



(264) 



(265) 



(266) 



Needless to say, all of these observables have to be evaluated in the absolute 
minimum of the effective potential with respect to m, i.e., by solving the gap 
equation. 

The phase diagram which one obtains in this way is depicted in Fig. p^ , 
adapted from Ref. 95. The line AB is a critical line of second order transitions 
(Voff changes from one minimum to a maximum and a minimum). The point 
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B is a tricritical point located at 1/(3 — 0.3183, /i — 0.6082 (all numbers here 
are in units of mo); it separates the second order line from a first order line BD 
along which Veff has two degenerate minima and a maximum. The endpoint 
D lies at fi — 1/V2 where the T = phase transition occurs, cf. Sec. 5.1 



Lines BC and BE are boundaries of metastability; when crossed, the potential 
acquires or loses a second minimum; point C is at = 1/2. In region OABD, 
chiral symmetry is broken and the quarks are massive; the outside region has 
unbroken chiral symmetry and massless fermions. 

In order to highlight the different behavior of the system when going 
through a first and second order phase transition, respectively, we have il- 
lustrated in Figs. ^ and ^ the evolution of VeS along the ^ = and T = 
axes (this corresponds to the special cases discussed in Sees. ^ and ^, respec- 
tively). Also shown is the behavior of the physical fermion mass. In this way, 
the unstable HF solution which we have found at finite density (with contin- 
uously decreasing mass, see Sec. 5.1) can be identified as corrsponding to a 
local maximum of Veff. In preparing Fig. 26, we have taken advantage of the 
fact that the integration in Eq. ( ^63| ) can be carried out at T = with the 
result 



}4ff 
N 



■In 



An 



(267) 



In 



6. 2 Chiral spiral at finite temperature 

We have seen before that the translationally non-invariant finite density solu- 
tion (constant baryon density, but oscillating condensates in the chiral limit) 
of the HF equation at T = can be obtained from the vacuum solution by a 
space dependmt chiral rotation of the quarks, cf. Sec. 3^. As already notec 
by Kao et aZJHJ for the massless Schwinger model and by Sadzikowski et alH^ 
for the a- and NJL models, the same procedure works at finite temperature as 
well. We show that the Fermi gas solution at zero chemical potential and finite 
temperature can be mapped onto the chiral spiral type solution at finite fi and 
the same T. This will allow us to follow the fate of the crystalline solution as 
we heat up baryonic matter. 

We start from the HF equation at finite T, but neither invoking trans- 
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Figure 25: Upper graph: Effective potential 
density as a function of dynamical fermion 
mass m at fi = for various temperatures 
below and above Tc ~ 0.567. Lower graph: 
m as a function of T (i.e., position of mini- 
mum in upper curves). Second order phase 
transition. 




Figure 26: Upper graph: Effective potential 
density as a function of dynamical fermion 
mass m at T = for various chemical po- 
tentials below and above fic = 1/ \/2. Lower 
graph: m as a function of fi (dashed line: 
position of maximum, corresponding to Fig. 
14; solid line: position of minimum in upper 
curves). First order phase transition. 



lational invariance nor assuming that (qij^q) vanishes, 

{^'t^ - ((g(2:)g(a;))7° + {qix)i^'q{x))i^')^ ^n{x) = c.„^„(.t) . (268) 

As before, temperature and chemical potential only show up in the self-consis- 
tency relations which now read 



{q{x)q{x)) = N^ipnix)fni3 



1 



n 

(q{x)i-f^qix)) = Nj2^nix)h'iPnix) , . ■ (269) 

n 
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Guided by the T = results (of. Sec. 5.3.1), we make the following ansatz for 
the a;-dependence of the scalar and pseudo-scalar condensates, 

Tfl r Tfl 

{qix)q{x)) = -cos{2fix) , {q{x)i-f'> q{x)) ^—sm{2^ix) . (270) 

9 9 

The HF equation for this particular potential can be turned into 

7'T^ + "^e'^"^'7°e-'''"'^'W„(x) =cu„^„(x) (271) 
1 ox J 

where we have used the fact that 



7° \ g-ipx^^ ^ f cos 2iix - sin 2^a; \ f 7° 
i'-f^ I \ sin2/ia; cos2/ia; J \ i-/^ 



(272) 



A local chiral rotation of the form 

^n{x) ^ e'''^''\n{x) (273) 

then eliminates the x-dependence from the HF potential, 

Id \ 
7^T7^+/i + "^7° Xn(a;) =w„Xn(a:) , (274) 



i dx 

and enables us to solve the HF equation trivially, 

Xn{x) = e''"'u{k) (w„ = e(fc) + /i) , 

Xn{x) = e'^'^vik) {iOn = -e{k) + fi) , (275) 



with e(fc) = \/W + rn^ and u{k),v{k) free, massive spinors. Since 

Xi{x)h'xn{x) = , (276) 
the two self-consistency conditions (|269|) collapse into the single condition 



^ -ExU^)x»(^) ^,(.„i,)^^ . (277) 



The additive term /i in cj„, Eq. ( |275| ), just cancels the chemical potential and 
we have mapped our problem exactly onto the /i = 0, translational invariant 
HF problem at finite T, 

~ ^ = I ' e„ = t^„-M = ±Vfc2 + m2 . (278) 

^ n 
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The parameter m therefore has to be identified with the physical fermion mass 
of the translationally invariant vacuum at /i = and the same finite T. We can 
thereby construct a translationally broken solution of the HF problem for any 
/i); it can be pictured as a chiral spiral whose radius shrinks with increasing 



temperature until it disappears at Tc = moe /tt, cf. Sec. 4.2 



Let us determine the free energy of this new phase. The density of the 



effective potential per color is still given by eq. (258), or equivalently by 



^ = ^ - — y In f 1 + e-^(-.-^)) . (279) 

The double counting term depends only on the radius of the chiral circle and 
is therefore unchanged (cf. Eq. (|270|)), 



(qq)' + {qij'qf = ^ . (280) 
9 

In view of the relation oj^ — n = £«, it appears superficially as if nothing 
depended on and we got the same result as for the translationally invariant 
system at zero chemical potential. As we know from the T — Q case, this 
argument is too rough — the /Lt-dependence is intimately related to the UV 
regularization. Let us impose an UV-cutoff on the single particle energies at 
the bottom of the Dirac sea, 

w„ > -A/2 ^ |fc|<A/2 + ^. (281) 

The entire /i-dependence of Vcff arises from this shift in the cutoff. We can 
then easily evaluate the difference between the free energy density of the "chiral 
spiral" at (/?, /J,) and the standard solution at (/3, 0), 



Voff(/3,M) 



N 



Veff(/3,0) 



N 

spir 



stand J-(A/2+/j) 27r 



A/2+A' 



dfc 



A/2 



J -A/2 



'-A/2 
(A + 



27r 



(282) 



This result does not depend on temperature, a feature characteristic for an UV 
effect. The fermion density per color can be obtained by differentiation with 
respect to /i; we find the same result as at zero temperature, 

P _ 9 Veff _ A /i 

iV-"5^^-2^ + ^' ^^^^^ 
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hence, with the definition of the Fermi momentum adopted in Sec. 



5.3.1 



^l^Pf (284) 
(the term ~ A in Eq. (283) is the fermion density of the Dirac sea). Summa- 



rizing, we can now compare the effective potentials for two solutions of the HF 
equation at finite (/3,/x). The Fermi gas solution with unbroken translational 



invariance has already been given in Eq. (263). The crystalline phase on the 
other hand yields 



N 





^m2\ 






An ^ 




spir 



4tt f3 J 2tt 



^ /'^'=ln(l + e-'''W)-^ , (285) 



where we have again dropped the term — A/i/(27r). In both cases, one still has to 
minimize the effective potential with respect to m. The depth of the respective 
minima (or physically minus the pressure) allows us to decide which phase is 
the stable one. The chiral spiral is favored everywhere over the homogeneous 
phase, below the critical temperature. This can be shown by a straightforward 
numerical evaluation, but also deduced analytically from the simple relation 

Vcff (/?, M)Lph. = Veff(/3, 0)-N^. (286) 



As the additional term in Eq. (286) does not depend on m, both Voff (/?, A')lspir 



and Veff(/3,0) are minimized by the same m 

mspir(/3,Ai) =m(/3,0) . (287) 

The free energy density of the spiral at P and /i evaluated at its minimum 
TOspir(/3, /i) is thus equal to 

Voff(/3,/i)lsph-,™=™.,„ = Veff(/3,0)|„(^^o) - . (288) 

'^csiP, ^J')\,n=m{t3 ^ monotonously decreasing function of /i at fixed /3 (its 
derivative with respect to fi at fixed /3 is just ~p). Furthermore from the phase 
diagram (Fig. |2^) it is clear that for any fixed P there is a critical chemical 
potential /Zc, where chiral symmetry gets restored. Above this critical chemical 
potential, Ves{P,^J,) evaluated at its minimum just behaves like the energy 
density of free fermions in 1+1 dimensions 

Veff(/3,M)U=™(/3,^) = -n!^ for/z > /X, . (289) 
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Figure 27: Phase diagram of chiral GN model, based on spatially inhomogeneous solution. 
The horizontal line at Tc ~ 0.567 is a second order line of chiral symmetry restoration. The 
dependence of radius and helix angle of the chiral spiral on /i and T is indicated qualitatively 
by means of the inserted drawings. 



Comparing Eq. ( P88|) to Veff (/?, Ai)lm=m(/3 fj.) clear that the chiral spiral is 
always lower in free energy than the homogeneous phase. 

Eq. ( 286 ) is also useful if one wants to evaluate the thermodynamic ob- 
servables for the crystal phase 

^ TT 

KAM)lsph. - K/3,0), 

£(/?,M)Lph. = e{P,0) + N^. (290) 

(At fi — 0, there is no distinction between the Fermi gas and the other phase.) 
Since the translationally broken phase is the thermodynamically stable phase 
wherever it exists, the phase diagram of the chiral GN model has only one 
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horizontal second order line at T — Tc where chiral symmetry is restored and 
the radius of the chiral spiral shrinks to zero (Fig. no phase transition 

occurs if we increase the chemical potential. Note that our results are rather 
different from a recent study of a non- uniform chiral phase in 3+1 dimensional 



effective field theories (cr-model, NJL modelled). In our language, these authors 
use the radius and the helix angle of the chiral spiral as independent variational 
parameters. Since in 1+1 dimensions baryon number is topologically connected 
to winding number, we do not seem to have this option but can vary only with 
respect to the radius (or m). 

As a last remark, we wish to comment on the original GN model with 
only discrete chiral symmetry. The chiral spiral does not exist here and one 
might think that our analysis does not have anything to say about it. However, 



the criticism of Sec. 5.1 also applies here. As the non-chiral GN model has 
(massive) bound baryons, the low density behavior of the energy obtained in 
standard HF approximation cannot be correct and the phase diagram has to 
be reconsidered. This is a harder task than for the chiral model since here 
dynamics is not completely dominated by chiral symmetry. It is quite possible 
that the first order phase transition survives, but this still has to be worked 
out in detail. 



7 Closing remarks 

In this paper, we have illustrated the kind of insights one gets from studies of 
exactly soluble field theories by means of two classic examples, the GN and the 
't Hooft model. The undeniable success of these models stems from a fortunate 
combination of two very bold substitutions: Replacing the number of spatial 
dimensions by 1 and the number of colors by oo (ironically, both are 3 in real 
QCD). Without the first choice, there is no solvability; without the 2nd, there 
are hardly enough interesting phenomena due to rigorous no-go theorems. If 
one is willing to swallow both of these bitter pills, one is awarded with a striking 
number of analogies between these toy models and the real world which have 
led to a long list of QCD-motivated applications. These analogies include such 
central features of QCD as asymptotic freedom, dimensional transmutation, 
chiral symmetry breaking, and confinement. In view of the fact that these 
models have the simplest Lagrangians of certain types which one can write 
down and are not contrived in any sense, this is quite amazing and underlines 
their pedagogical value for QCD related questions. Nowadays, it may well be 
that the study of lower dimensional field theories is a better training ground 
for students interested in QCD than a standard field theory course based on 
Feynman diagrams and perturbative QED. 
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Throughout this work we have tried to emphasize the fact that the solution 

of these models can be found by rather elementary means we have delib- 
erately chosen the relativistic HF approach as our main tool. This does not 
mean though that there are no subtleties: We only mention the emergence of 
baryon number from the UV regularization, or how zero mode gluons affect the 
thermodynamics of quarks. Such features make the use of analytical methods 
particularly important and set the scale for the even more subtle properties to 
be expected from 4-dimensional QCD. 

As far as our more specific subject — finite temperature and density 
behavior — is concerned, we have found that even after more than 25 years 
these models are still good for some surprises. The strongest impression for us 
is the way in which chiral symmetry governs not only the hadron spectrum, 
but also the phase structure at finite density and temperature. Confinement 
is most clearly seen in the temperature dependence of observables, whereas 
the dependence on the chemical potential is dominated by chiral symmetry. 
Other non-trivial insights are related to the realization of the Skyrme picture 
in a precise fashion, a triple point in the phase diagram of the non-chiral 
GN model and inhomogeneous chiral condensates. As a by-product, we have 
resolved some puzzles in the literature and identified issues where further work 
is needed, for instance the phase diagram of the non-chiral GN model. 

In conclusion, we see no indication that models such as the GN or 't Hooft 
model are anywhere close to having exhausted their potential as a theoretical 
laboratory for QCD. 

Acknowledgements 

We are grateful to Misha Shifman for the invitation to contribute to this 
Festschrift in honor of Boris loffe, thereby forcing us to clarify some issues 
before a certain deadline. M.T. thanks his collaborators, colleagues and stu- 
dents for most enjoyable and useful interactions on 2d field theory, in particular 
M. Burkardt, V. Dolman, M. Engelhardt, O. Jahn, M. Kurz, F. Lenz, S. Levit, 
R. Oderkerk, K. Ohta, R. Pausch, C. Satzinger, B. Schreiber, M. Seeger, A. 
Smilga, D. Stoll, R. Weingartner, S. Worlen, K. Yazaki, and I. Zahed. 

References 

1. N. Seiberg and E. Witten, Nucl. Phys. B 426, 19 (1994). 

2. J. Schwinger, Phys. Rev. 128, 2425 (1962). 

3. D.J. Gross and A. Neveu, Phys. Rev. D 10, 3235 (1974). 

4. G. 't Hooft, Nucl. Phys. B 75, 461 (1974). 

5. D.R. Hartree, Proc. Camb. Phil. Soc. 24, 89 (1928). 



85 



6. V.A. Fock, Z. Phys. 61, 126 (1930). 

7. M. Blau and G. Thompson, Trieste HEP CosmoL, 0175-244 (1993). 

8. E. Abdalla and M.C.B. Abdalla, Phys. Rep. 265, 253 (1996). 

9. S.J. Brodsky, H.-C. Pauli, and S.S. Pinsky, Phys. Rep. 301, 299 (1998). 

10. N. Nakanishi and I. Ojinia, Covariant operator formalism of gauge theo- 
ries and quantum gravity, World Scientific 1991. 

11. S. Coleman, Comm. Math. Phys. 31, 259 (1973). 

12. C. Itzykson and J.-B. Zuber, Quantum field theory, McGraw-Hill 1980. 

13. N.D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 (1966). 

14. K. Huang, Quantum Field Theory. From Operators to Path Integrals, 
John Wiley 1998, ch. 18. 

15. E. Witten, Nucl. Phys. B 145, 110 (1978). 

16. I. Affleck, Nucl. Phys. B 265, 448 (1986). 

17. J.M. Kosterlitz and D.J. Thouless, J. Phys. C 6, 1181 (1973). 

18. J.M. Kosterlitz, J. Phys. C 7, 1046 (1974). 

19. V.L. Berezinsky Sov. Phys. JETP 32, 493 (1971). 

20. N.S. Manton, Ann. Phys. (N.Y.) 159, 220 (1985). 

21. S.G. Rajccv, Phys. Lett. B 212, 203 (1988). 

22. J.E. Hctrick and Y. Hosotani, Phys. Lett. B 230, 88 (1989). 

23. J.E. Hctrick, Int. J. Mod. Phys. A 9, 3153 (1994). 

24. F. Lenz, H.W.L. Naus, and M. Thies, Ann. Phys. (N. Y.) 233, 51 (1994). 

25. J.A. Minahan and A.P. Polychronakos, Phys. Lett. B 326, 288 (1994). 

26. B. Svetitsky Phys. Rep. 132, 1 (1986). 

27. F. Lcnz, M. Shifman, and M. Thics, Phys. Rev. D 51, 7060 (1995). 

28. M. Engclhardt and B. Schrcibcr, Z. Phys. A 351, 71 (1995). 

29. A. Dhar, G. MandaL and S.R. Wadia, Nucl. Phys. B 436, 487 (1995). 

30. M. Sccgcr and M. Thics, Phys. Rev. D 58, 027701 (1998). 

31. A. Dhar, G. Mandal, and S.R. Wadia, Phys. Lett. B 329, 15 (1994). 

32. A. Dhar, P. Lakdawala, G. Mandal, and S.R. Wadia, Int. J. Mod. Phys. 
A 10, 2189 (1995). 

33. W.E. Thirring, Ann. Phys. (N.Y.) 3, 91 (1958). 

34. G. 't Hooft, Nucl. Phys. B 72, 461 (1974). 

35. S. Coleman, Erice lecture 1979, in: Aspects of symmetry, Cambridge 
University press 1985, p. 351. 

36. E. Witten, Nucl. Phys. B 160, 57 (1979). 

37. A.V. Manohar, in: Probing the Standard Model of Particle Interactions, 
F. David and R. Gupta, eds. (Les Houches 1997). 

38. Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1960); 124, 246 
(1961). 

39. A.B. Zamolodchikov and A.B. Zamolodchikov, Phys. Lett. B 72, 481 



86 



(1978). 

40. N. Andrei and J.H. Lowenstein, Phys. Rev. Lett. 43, 1698 (1979). 

41. N. Andrei and J.H. Lowenstein, Phys. Lett. B 90, 106 (1980); B 91, 401 
(1980). 

42. R.F. Dashen, B. Hasslacher, and A. Neveu, Phys. Rev. D 12, 2443 
(1975). 

43. J. Feinberg, Phys. Rev. D 51, 4503 (1995). 

44. L.L. Salcedo, S. Levit, and J.W. Negele, Nucl. Phys. B 361, 585 (1991). 

45. F. Lenz, M. Thies, S. Lcvit, and K. Yazaki, Ann. Phys. (N. Y.) 208, 1 
(1991). 

46. R. Pausch, M. Thies, and V.L. Dolman, Z. Phys. A 338, 441 (1991). 

47. D.J. Gross and F. Wilczek, Phys. Rev. Lett. 30, 1343 (1973). 

48. H.D. Politzer, Phys. Rev. Lett. 30, 1346 (1973). 

49. W. Wetzel, Phys. Lett. B 153, 297 (1985). 

50. C. Luperini and P. Rossi, Ann. Phys. (N.Y.) 212, 371 (1991). 

51. S. Shei, Phys. Rev. D 14, 535 (1976). 

52. J. Feinberg and A. Zee, Phys. Lett. B 411, 134 (1997). 

53. J. Feinberg and A. Zee, Phys. Rev. D 56, 5050 (1997). 

54. A.C. Scott, F.Y.F. Chu, and D.W. McLaughlin, Proc. of the LEEE 61, 
1443 (1973). 

55. C.G. Callan, S. Coleman, D.J. Gross, and A. Zee, unpublished (referred 
to by Ref. 42). 

56. M. Thies and K. Ohta, Phys. Rev. D 48, 5883 (1993). 

57. M.B. Einhorn, Phys. Rev. D 14, 3451 (1976). 

58. C.G. Callan, N. Coote, and D.J. Gross, Phys. Rev. D 13, 1649 (1976). 

59. M.B. Einhorn, S. Nussinov, and E. Rabinovici, Phys. Rev. D 15, 2282 
(1977). 

60. J. Ellis, Acta Phys. Polon. B 8, 1019 (1977). 

61. J.H. Weis, Acta Phys. Polon. B 9, 1051 (1978). 

62. T.T. Wu, Phys. Lett. B 71, 142 (1977). 

63. I. Bars and M.B. Green, Phys. Rev. D 17, 537 (1978). 

64. M. Li, Phys. Rev. D 34, 3888 (1986). 

65. M. Li, L. Wilets, and M.C. Birsc, J. Phys. G 13, 915 (1987). 

66. S. Dalley and L Klebanov, Phys. Rev. D 47, 2517 (1993). 

67. G. Bhanot, K. Demeterfi, and L Klebanov, Phys. Rev. D 48, 4980 
(1993). 

68. I. Bigi, M. Shifman, N. Uraltsev, and A. Vainshtein, Phys. Rev. D 59, 
054011 (1999). 

69. B. Grinstein and R.F. Lebed, Phys. Rev. D 59 , 054022 (1999). 

70. M. Burkardt and N. Urahsev, ^iep-ph/0005278 . 



87 



71. A.R. Zhitnitsky, Phys. Lett. B 165, 405 (1985). 

72. M. Burkardt, Phys. Rev. D 53, 933 (1996). 

73. T.H. Koopmans, Physic a 1, 104 (1933). 

74. V. Schon and M. Thies, |hep-th/0003T95[ 



75. A.L. Fetter and J.D. Walecka, Quantum theory of many-particle systems, 
McGraw-Hill 1971. 

76. M. Burkardt, Multiquarksysteme in der 1+1 dimensionalen QCD, phd 
thesis, Universitat Erlangcn-Niirnberg 1989. 

77. M. Durgut, Nucl. Phys. B 116, 233 (1976). 

78. D. Amati and E. Rabinovici, Phys. Lett. B 101, 407 (1981). 

79. D. Amati, K.-C. Chou, and S. Yankielowicz, Phys. Lett. B 110, 309 
(1982). 

80. R.P. Oderkerk, Relativistic hadrons in an external field, phd thesis. Free 
University of Amsterdam (VU) 1992. 

81. P. Ring and P. Schuck, The nuclear many-body problem, Springer, New 
York 1980. 

82. C. Satzinger, Wiederherstellung der chiralen Symmetric in QCD2, Diplo- 
marbeit, Universitat Erlangen-Niirnberg 1991. 

83. P.H. Damgaard, H.B. Nielsen, and R. Sollacher, Nucl. Phys. B 385, 227 
(1992). 

84. T.H.R. Skyrme, Proc. Roy. Soc. Lond. A 260, 127 (1961). 

85. M. Gell-Mann, R.J. Oakes, and B. Renner, Phys. Rev. 175, 2195 (1968). 

86. F. Lenz and M. Thies, Ann. Phys. (N.Y.) 268, 308 (1998). 

87. D.J. Toms, Phys. Rev. D 21, 928 (1980). 

88. V. Koch, E.V. Shuryak, G.E. Brown, and A.D. Jackson, Phys. Rev. D 
46, 3169 (1992). 

89. A.P. Young, Nucl. Phys. B 42 (Proc. SuppL), 201 (1995). 

90. S.L. Sondhi, S.M. Girvin, J.P. Carini, and D. Shahar, Rev. Mod. Phys. 
69, 317 (1997). 

91. J.I. Kapusta, Finite-temperature field theory, Cambridge Univesity Press 
1989. 

92. A. Barducci, R. Casalbuoni, M. Modugno, G. Pettini, and R. Gatto, 
Phys. Rev. D 51, 3042 (1995). 

93. R.F. Dashen, S.K. Ma, and R. Rajaraman, Phys. Rev. D 11, 1499 
(1975). 

94. F. Karsch, J. Kogut, and H.W. Wyld, Nucl. Phys. B 280, 289 (1987). 

95. U. Wolff, Phys. Lett. B 157, 303 (1985). 

96. T. Eguchi and H. Kawai, Phys. Rev. Lett. 48, 1063 (1982). 

97. B.J. Harrington and A. Yildiz, Phys. Rev. D 11, 779 (1975). 

98. L.D. McLerran and A. Sen, Phys. Rev. D 32, 2794 (1985). 



88 



99. T.H. Hansson and I. Zahed, Phys. Lett. B 309, 385 (1993). 

100. M. Engclhardt, Phys. Lett. B 355, 507 (1995). 

101. M. Thies, Nucl. Phys. A 546, 233c (1992). 

102. V. Schon and M. Thies, Phys. Lett. B 481, 299 (2000). 

103. J. Gasser and H. Leutwyler, Phys. Lett. B 184, 83 (1987). 

104. P. Gerber and H. Leutwyler, Nucl. Phys. B 321, 387 (1989). 

105. P. Zhuang, J. Hiifner, and S.P. Klevansky, Nucl. Phys. A 576, 525 
(1994). 

106. S.P. Klevansky, Rev. Mod. Phys. 64, 649 (1992). 

107. A. Barducci, R. Casalbuoni, M. Modugno, G. Pettini, and R. Gatto, 
Mod. Phys. Lett. A 11, 1579 (1996). 

108. A. Barducci, M. Modugno, G. Pettini, R. Casalbuoni, and R. Gatto, 
Phys. Rev. D 59, 114024 (1999). 

109. A. Barducci, R. Casalbuoni, M. Modugno, G. Pettini, and R. Gatto, 
Phys. Rev. D 55, 2247 (1997). 

110. T.F. Treml, Phys. Rev. D 39, 679 (1989). 

111. A. Chodos, R.L. Jaffe, K. Johnson, C.B. Thorn, and V.F. Weisskopf, 
Phys. Rev. D 9, 3471 (1974). 

112. M. Alford, K. Rajagopal, and F. Wilczek, Phys. Lett. B 422, 247 

(1998). 

113. I. Klcbanov, Nucl. Phys. B 262, 133 (1985). 

114. W.L. McMillan, Phys. Rev. B 16, 4655 (1977). 

115. G. Thcodorou and T.M. Rice, Phys. Rev. B 18, 2840 (1978). 

116. K. Takayama and M. Oka, Nucl. Phys. A 551, 637 (1993). 

117. Handbook of Mathematical Functions, M. Abramowitz and I. A. Stegun, 
eds., Dover, Now York 1970. 

118. W. Fischlcr, J. Kogut, and L. Susskind, Phys. Rev. D 19, 1188 (1979). 

119. E. Wigncr, Phys. Rev. 46, 1002 (1934). 

120. H.J. Schulz, Phys. Rev. Lett. 71, 1864 (1993). 

121. Y.-C. Kao and Y.-W. Lcc, Phys. Rev. D 50, 1165 (1994). 

122. H.R. Christiansen and F.A. Schaposnik, Phys. Rev. D 53, 3260 (1996). 

123. H.R. Christiansen and F.A. Schaposnik, Phys. Rev. D 55, 4920 (1997). 

124. M. Kutschcra, W. Broniowski, and A. Kotlorz, Nucl. Phys. A 516, 566 
(1990). 

125. D.V. Doryagin, D.Yu. Grigoriev, and V.A. Rubakov, /ni. J. Mod. Phys. 
A 7, 659 (1992). 

126. A.W. Overhauscr, Adv. in Phys. 27, 343 (1978). 

127. E. Shuster and D.T. Son, Nucl. Phys. B 573, 434 (2000). 

128. B.-Y. Park, M. Rho, A. Wirzba, and I. Zahed, Phys. Rev. D 62, 034015 
(2000). 



89 



129. A. Chodos and H. Minakata, Phys. Lett. A 191, 39 ( 1994). 

130. M. Sadzikowski and W. Broniowski, |hcp-ph/0003282[ 



90 



